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motion  and  that  of  the  gun  are  independent,  and  they  may  be  analyzed  by  the 
dynamical  theory  of  a  single  rigid  body.  The  motion  of  the  projectile  with 
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respect  to  the  tube  is  regarded  as  known.  Also,  the  base  pressures  on  the 
projectile  and  the  tube  are  regarded  as  known  functions  of  time.  Geometrically, 
the  projectile  is  regarded  as  a  body  of  revolution  that  fits  snuggly  in  the  tube, 
so  that  there  is  no  balloting.  However,  a  lopsided  distribution  of  mass  density 
in  the  projectile  is  admitted  in  Sections  1,  2,  and  3.  In  Section  4,  only  a 
balanced  projectile  is  considered.  In  applications  of  the  momentum  principles, 
the  momentum  of  gases  in  the  tube  is  neglected.  It  can  be  introduced  empirically 
by  augmenting  the  mass  of  the  projectile  by  a  portion  of  the  mass  of  the  charge. 
This  concept  of  effective  mass  has  a  precedent  in  other  branches  of  fluid 
mechanics. 

Section  1  treats  a  gun  that  is  completely  unsupported.  The  free  gun  is  an 
extreme  case;  since  the  recoil  mechanism  and  the  trunnion  would  restrain  the 
motion.  Section  2  treats  another  extreme  case;  namely,  a  rigid  gun  that  is 
completely  fixed.  In  this  case,  there  is  no  problem  of  motion,  since  the  motion 
of  the  projectile  is  presumed  to  be  given.  Consequently,  the  analysis  deals  only 
with  the  forces  and  moments  of  interaction  between  the  gun  and  the  projectile. 
Section  3  treats  the  motion  of  a  gun  t1  at  translates  freely  along  a  guide,  but 
it  is  constrained  against  rotation.  Section  4  treats  a  gun  with  a  trunnion  and 
a  recoil  mechanism  that  reacts  with  a  force  that  is  an  arbitrary  function  of  the 
displacement  and  velocity  of  recoil.  The  motion  is  constrained  to  take  place  in 
a  plane. 

Vectorial  mechanics  is  used  in  Sections  1,  2,  and  3.  In  Section  4,  the 
Lagrange  equations  (scalar  form)  are  used. 
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SECTION  1 

FREE  RIGID  GUN  TUBE  WITH  ECCENTRIC  BREECH  AND 
DYNAMICALLY  UNBALANCED  PROJECTILE 

1.1  INTRODUCTION 

The  motion  of  a  free  rigid  gun  in  a  gravitationless  environment 
is  investigated  in  this  section.  Also,  the  force  and  the  moment  that  the 
projectile  exerts  on  the  gun  are  derived.  Gravity,  acting  on  a  free  gun, 
would  cause  the  center  of  mass  of  the  system  to  descend  with  constant 
acceleration  g.  The  supports  of  an  actual  gun  prevent  a  free  fall  of  the 
system.  Consequently,  it  is  questionable  whether  inclusion  of  gravity  in 
the  analysis  would  make  the  mathematical  model  more  realistic.  Like  gravity, 
the  resistance  of  air  in  the  tube  ahead  of  the  projectile  is  an  external 
force.  Its  magnitude  can  be  calculated  by  gas  dynamics.*  However,  in  this 
chapter,  the  air  resistance  on  the  projectile  is  disregarded.  The  system, 
as  it  is  conceived,  accordingly  has  no  external  forces  acting  on  it.  Since 
the  relative  motion  of  the  projectile  with  respect  to  the  tube  is  assumed 
to  be  given,  the  motion  of  the  system  is  determined  by  the  laws  of  conserva¬ 
tion  of  rectilinear  and  angular  momentum. 

1.2  TERMINOLOGY  AND  NOTATIONS 

The  gun  tube  and  the  breech  block  together  form  a  rigid  body 
called  the  "gun.”  The  system  consists  of  the  gun  and  the  projectile  in  the 
tube.  The  point  on  the  geometric  axis  of  the  projectile  that  lies  closest 
to  the  center  of  mass  of  the  projectile  is  called  the  "geometric  center"  of 
the  projectile.  The  location  of  the  center  of  mass  of  the  projectile  at  the 
initial  instant  is  called  the  "starting  point  of  the  center  of  mass  of  the 
projectile." 

A  bar  over  a  letter  denotes  a  vector.  A  caret  over  a  letter 
denotes  a  unit  vector.  A  dot  over  a  letter  denotes  the  derivative  with 
respect  to  time  t.  Ignition  occurs  at  the  instant  t  =  0.  The  following 
notations  are  illustrated  by  Figure  1. 

Point  0  is  the  center  of  mass  of  the  whole  system. 

Point  P  is  the  center  of  mass  of  the  gun. 

R.  Courant  and  K.  Friedrichs ,  "Supersonic  Flow  and  Shock  Waves  (U) 

Chap.  LTI ,  Interscience  Publishers ,  New  York,  1948. 
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Point.  Q  is  the  center  of  mass  of  the  projectile  at  time  t. 

Point  Q'  is  the  geometric,  center  of  the  projectile  at  time  t. 

Point;  Qq  is  the  starting  point,  of  the  center  of  mass  of  the 

project  tie. 

Point  Qy ’  is  the  initial  location  of  point  Q' . 

f*  denotes  a  Galilean  reference  frame;  e.g.,  the  earth. 

x , y , are  rectangular  coordinates  fixed  in  frame  F. 

x'/jVic'  are  unit  vectors  along  the  axes  x,y,z. 

A  A  A 

i,j,k  are  unit  vectors  along  the  principal  axes  of  inertia  of  the 
gun  through  point  P. 

a,6,c  are  unit  vector's  along  the  principal  axes  of  inertia  of  the 
projectile  through  point  Q.  None  of  them  is  necessarily  parallel  to  the 
axis  of  the  tube. 

AAA 

a0’k(),c0  arc  tmit  vectors  along  the  principal  axes  of  inertia  of 
the  projectile  through  point  Q0,  when  the  projectile  is  in  the  starting 
position. 

v  is  a  unit  vector  along  the  axis  of  the  tube. 

R()  is  the  vector  oi\ 
is  the  vector  oq. 

r  is  the  vector  PQ;  r  =  Uj  -  R  . 

e  is  the  vector  I’Q^ 1 . 

e  is  the  vector  Q^Q, 
is  the  vector 

s  is  the  distance  that  the  projectile  has  traveled  relative  to  the 
tube  at  time  t. 


X  is  the  angle  through  which  the  projectile  has  turned  relative 
to  the  tube  at  time  t. 

V 

w  r-  x  is  the  spin  (angular  velocity)  of  the  projectile  relative 
to  the  tube  at  time  t.. 

M  is  the  mass  of  the  gun. 

m  is  the  mass  of  the  projectile. 

*1’*2’*3  Urc  Pr i ncipal  moments  of  inertia  of  tho  gun  with 
respect  to  its  center  of  mass  P. 

ij.i^.i-  are  the  principal  moments  of  inertia  of  the  projectile 
with  respect  to  its  center  of  mass  Q. 
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W  is  the  angular  velocity  of  the  gun  relative  to  frame  /•. 
w  is  the  angular  velocity  of  the  projectile  relative  to  frame  t: . 
WrW2,W_  are  components  of  W,  defined  by  W  -  iW^  jW.,  +  kW^. 

wl,w2>w3  are  comPonents  o1'  w,  defined  by  w  -■  aw ^  +  bw^  +  cw „ . 

rl*r’,r3  urc  coml,onents  of  r,  defined  by  r  =  ii^  +  ji'2  +  Icr^. 

Xp,Y^,2p  arc  the  (x,y,z)  coordinates  of  the  center  of  mass  P 

of  the  gun;  i.c.,  R0  =  i*X0  +  ]’Y0  +  k*Z(). 

el,e2,63  aro  comP°non!ts  the  vector  e",  defined  by  e  =  ie^  + 

K  *  k«,. 

ot,(3,y  are  direction  cosines  of  the  axis  of  the  tube,  defined  by 

A  A  A  A 

v  =  ia  +  j(3  +  ky. 

i||  _  A 

ej,e2,Gj  are  components  of  vector  defined  by  =  ie^  + 

je2  +  key 

A  A  A 

i,mi,ni  are  direction  cosines  of  vectors  i,  j,  k  with  respect 
to  axes  (x,y,2),  defined  by  the  matrix 

AAA 

i'  j*  k' 

i  i»j 

3  l2  ", 

k 


ai,bi,ci  are  direction  cosines  of  vectors  (a,b,c),  defined  by  the 


matrix 


AAA 

i  j  k 

u  ax  a2  a3 

b  b,  b2  b3 


c  Cj  c2  c3 


^  A  A 

ot.  ,y.  are  direction  cosines  of  vectors  (a  ,1)  ,c  ),  defined  by 
1X1  J  u  u 


the  matrix 


A  A  A 


i  j  k 
\)  al  B1  Y1 

b0  a2  Y2 
a0  U3  6  3  Y3 
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iT  is  the  angular  momentum  of  the  gun  about  point:  0. 
h  is  the  angular  momentum  of  the  projectile  about  point  0. 

Py(t)  is  the  pressure  of  the  gases  on  the  breoch. 

Pj(t)  is  the  pressure  of  the  gases  on  the  base  of  the  projectile. 
A  is  the  cross-sectional  area  of  the  bore. 

F  is  the  resultant  force  exerted  by  the  projectile  on  the  tube. 

M  is  the  moment  about  the  center  of  mass  of  the  gun  of  the 
contact  forces  that  the  projectile  exerts  on  the  tube. 

M'  is  the  moment  about  the  center  of  muss  of  the  projectile  of 
the  contact  forces  that  the  tube  exerts  on  the  projectile, 
g  is  the  scalar  acceleration  of  gravity, 
g  is  the  vector  acceleration  of  gravity  (directed  vertically 


downward) 


R  is  the  resisting  force  of  the  air  ahead  of  the  projectile 


"  =  M  +'  I 

arc  components  of  vector  M,  defined  by  M  =  rMj  +  jM^  + 

kM  y 

I'l  *  I  2  *  F-s  iU'°  c0inlX)n0llts  °f  vector  F,  defined  by  V  -  ifj  +  jl'2  + 
l<F 

3* 

f  =  if  +  jf  +  Kf  is  the  correction  to  force  F  to  account  for 

1  M  J 

unbalance  of  the  projectile. 

1.3  CONSlillVATTON  Of  MOMENTUM 

Since  there  are  no  external  forces,  the  center  of  mass  0  of  the 
system  remains  fixed  in  Frame  F.  Consequently,  in  this  section,  0  is  taken 
to  be  the  origin  of  coordinates  (x,y ,/.) .  At  time  t  ,  the  center  of  mass  Q  of 
the  projectile  lies  at  point  +  v  (figure  1).  Since  the  origin  0  is 

the  center  of  mass  of  the  system,  and  since  momentum  of  the  gases  is 
neglected , 


MRq  +•  mRj  =  0 


0.1) 


Consequently, 


MRQ  +  m(R()  +  r)  =  0 
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Therefore, 


m 


1 


m 


n  _  -mr  .  r  _  Mr 
0  “  M  +  m  *  1  "  M  +  m 


-  —  R 
m  K0 


(1.2) 


Equation  (1.2)  signifies  that  the  vectors  R^  and  are  collinear, 
Consequently,  point  0  lies  on  the  line  PQ.  Equation  (1.2)  yields 


_  jl  _  a.  Mm  — 

MRn  x  R.  +  mR.  x  R,  =  - r  x  r 

0  0  1  1  M  +  m 


(1.3) 


A  /. 


Equation  (1.3)  will  be  used  later. 

The  direction  cosines  of  vectors  (i,j,k)  with  respect  to  axes 
(x,y,z)  are  (JL,nt,n^)  (see  Notations).  Also 


«o  ■  i,x0  *  ^'Yo  *  k'zo 


Equation  (1.2)  yields 


R0  ■  *  *2  *  kr3> 


Also, 


A  A  A  A 

i  =  £^i  1  +  j  ’  +  n,  k 1 


£?i ’  +  m2j  *  +  n2k' 


A  A 


k  =  Jl^i '  +  m^j  •  +  n^k ' 


Consequently, 


X0  ”  M  +  m(rlAl  +  r2^2  +  r3£3) 
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ANGULAR  VELOCITY  OF  THE  PROJECTILE 

The  absolute  angular  velocity  of  the  projectile  is 

w  =  aw^  +  b«2  +  cWj 


Let  (a'.b'.c1)  be  an  orthogonal  triad  of  unit  vectors,  such  that  c’  =  v, 
where  v  is  the  unit  vector  along  the  axis  of  the  tube.  Then  w  may  be  written 
alternatively  as  follows: 

w  =  a'w^'  +  +  c'Wj* 

The  angular  velocity  of  the  projectile  relative  to  the  tube  is  a 
vector  of  magnitude  to  =  x  that  is  coaxial  with  the  tube.  The  variable  w  is 
the  spin;  it  is  regarded  as  a  known  function  of  t.  Accordingly,  the  cf 
component  of  the  absolute  angular  velocity  of  the  projectile  is 


w3*  =  w  +  W  *  c* 


where  W  is  the  angular  velocity  of  the  gun.  The  a'  and  b*  components  of  the 
absolute  angular  velocity  of  the  projectile  are  the  same  as  those  of  the  tube, 
since  balloting  is  excluded.  Accordingly, 


Wj'  =  a '  •  W  ;  w^*  =  S’  •  W 
* 

Consequently  , 

w  =  a’a’  •W  +  b’b'  •  W  +  c'c'  •  W  ♦  c'w 

The  notation  a'a'  •  W,  eta.,  is  short  for  a '(a1  •  W).  Since  a '  •  W  is  a 
scalar,  a'a’  •  W  denotes  multiplication  of  vector  a'  by  the  scalar  a'  •  W. 


However,  the  following  relationship  is  an  identity: 


a’a'  •  W  +  b'b'  •  W  +  c'c'  •  W  =  W 


A  -A 


Also,  c'  =  v.  Consequently, 


w  =  W  +  vti) 


Cl.  5 


Since, 


w  =  aw^  +  bw^  +  cw^  ;w1=w*a;w2=w*b;w3=w»c 


Also, 


A  A 


a  =  iaj  +  ja2  +  ka3 


A  A 


b  =  ibj  +  jb2  +  kb3 


A  A 


c  =  iCj  +  jc2  +  kc3 


Consequently,  Eq.  (1.5)  yields 


”l  "  Vl  *  a2W2  *  a3W3  *  (aal  *  6a2  *  Ya3)w 


w2  =  b1W1  +  b2W2  ♦  b3W3  +  (ab1  +  gb2  +  yb3)0) 


w.  =  CjWj  +  c2W2  ♦  c3W3  +  (ocj  +  6c2  *  yc3)a, 


(1.6 


1.5  THE  VECTORS  a,b,c,e 

AAA 

So  far,  the  vectors  a,  b,  c,  whose  componer  ppear  in  Eq.  (1.6) 
are  undetermined.  The  displacement  vector  field  of  a  ,  igid  body  that  under¬ 
goes  a  rotation  about  an  axis  that  is  oblique  to  the  coordinate  ax"  »  is 

AAA 

derived  in  Appendix  A.  Vectors  a,  b,  c  are  imbedded  in  the  projectile  and 


move  with  it.  Accordingly,  Eq.  (A-2)  in  Appendix  A  applies  to  them.  The 
vector  (T  in  Eq.  (A-2)  may  be  given  various  special  designations.  If  p^= 
a0>  a  =  R.  If  p  =  bQ,  b  =  R,  and  likewise  for  cQ  and  c.  Consequently  , 


A  A  A 


a  =  vxaQsinx  +  vv 


A  /V  A 


b  =  vxbQsinx  +  v'» 


/\  A  A 


c  =  vxcQsmx  +  w 


aQ(l  -  cosy)  +  aQcosx 
bQ(l  -  cosy)  +  bQcosx 
cQ(l  -  cosy)  +  c0cosx 


(1.7) 


It  can  be  shown  directly  from  Eq.  (1.7)  that  (a,b,c)  are  an  orthogonal  triad 
of  unit  vectors,  as  they  should  be.  Equation  (1.7)  also  yields 


A  A  A  A  A  A  A  A  A  A  A  A 

v*a=v»a0;v*b=v*b0;v*c=v*c0 


(1.8) 


Equation  (1.8)  reflects  the  fact  that  a  rotation  of  a  body  about  an  axis  does 
not  change  the  angle  between  the  axis  and  any  line  that  is  scribed  in  the 
body. 

The  direction  cosines  of  the  vectors  (a0,bp,c0)  are  given  by  the 
following  matrix: 


(1.9) 


The  direction  cosines  in  Eq.  (1.9)  are  known  constants.  Equations  (1.7)  and 
(1.9)  yield 


★ 

See  footnote  on  page 


13  and  Appendix  C. 


a  =  i  [C By x  -  YB^siny  +  aCcwij  +  $$1  +  yy^O  -  cosy)  +  o^cosy] 


i 

! 


s. 


I 

-^1 


*  j  [(ya1  -  ay^siny  +  Bfaoij  +  BBj  +  yy^Cl  -  cosy)  +  B^osy] 

+  k  [(aBj  -  Ba^siny  +  y(cta1  +  0B1  +  YYjMI  -  cosx)  +  y^osy] 

b  =  i  [(By2  -  yB2)sinx  +  a(aa2  +  BB2  +  yy2)(l  -  cosy)  +  o^cosy] 

+  j  [(Y“2  "  ay2)sinx  +  ^(aa2  +  3B2  +  yy2)(i  -  cosy)  +  B2cosx] 

+  £  [(aB2  -  Ba^siny  +  y(aa2  +  BB2  +  yy2)(l  -  cosy)  +  Y2cosy] 

c  =  i  [CBy3  -  yP3)siny  +  a(aa3  +  BB3  +  YY3)(1  -  cosy)  +  c^cosy] 

+  J  Ky^j  “  ay3)siny  +  3(aa3  +  B$3  +  YY3)(1  -  cosy)  +  BjCOsy] 

+  k[(ap3  -  Ba3)siny  +  y(a«3  +  $P3  +  yy3)(l  -  cosy)  +  y3cosy]  (1.10) 

Equation  (1.10)  determines  the  direction  cosines  (a.,b.,c.)  of  vectors 

AAA  1X1 

(a,b,c)  as  functions  of  t. 

If  p  is  set  equal  to  e  in  Eq.  (A-2) ,  R  =  e .  Since  e_  is 

A  ^  ^ 

perpendicular  to  v,  \>  •  £^=0,  Consequently ,  by  Eq.  (A-2), 

_  a  . ^ _ 

e  =  vxeQsiny  +  e0cosy  (1.11) 

Since  •  v  =  0,  it  follows  from  Eq.  (1.11)  that  e  •  v  -  0,  as  it  should  be. 

Also,  Eq.  (1.11)  yields 


where  e  is  the  eccentricity  of  the  projectile. 


1.6  THE  VECTOR  PRODUCT  r  x  r 

* 

The  vector  product  r  x  r  occurs  in  the  expression  for  the  angular 
momentum  of  the  system  about  point  0.  By  Figure  1, 

t  =  e  +  vs  +  e  (1.13) 

Hence , 

•  •  ^  ^  •  » 

r*  =  e  ♦  vs  +  vs  +  e 

_  A 

In  the  reference  frame  of  the  gun,  e  and  v  are  constant  vectors.  Conse¬ 
quently, 

.L  _  _  X  _ 

e=Wxe;v=Wxv  (1.14) 

Therefore, 

—  _  A 

r=Wxe+Wxvs+vs+e  (1.15) 

•  _  _ 

Since  e0  =  W  x  cQ,  Eq.  (1.11)  yields  (see  footnote,  page  13) 

i  A_  _  _A_  /N_  __  _ 

e  =  vW  •  e^siny  -  e^v  •  Wsinx  +  ojv  x  eQcosx  +  W  x  e^cosx  -  aie^siny  (1.16) 
Equations  (1.11),  (1.13),  (1.15),  and  (1.16)  yield 

_  A  A  _  _ 

r  =  e  +  vs  +  v  x  eQsinx  +  eQcosx  (1.17) 

r  =  W  x  e  +  W  x  vs  +  vs  +  vW  •  e^sinx  -  CqV  *  Wsinx 

A  _  _  _ _  _ 

+  wv  x  EgCosx  +  W  x  F-Ocosx  -  oje0sinx  (1.18) 

Equation  (1.18)  may  be  written  as  follows: 

4.  _  _  ^  #  -A—  _ 

r  =  W  x  r  +  vs  +  tov  x  e^cosy  -  we^sinx 


(1.19) 


Hence  (see  footnote,  page  13), 

r  x  r  =  Wr2  -  r?*W  +  rxvs  +  wvr  •  eQcosx  -  we0r  •  vcosx 


-  0) r  x  e^sinx 


(1.20) 


With 

r  =  ir^  +  jr2  +  kr3  ,  W  =  iW^  +  jW^  +  kWj  , 

0  =  ia  +  jB  +  kY  ,  and  eQ  =  iex  +  je2  +  ke3 


Eq.  (1.20)  yields 

7x7=  i{(r22  .  r/jKj  -  rjr2W2  -  <■  (yx2  -  BXj^ 

+  [a(r2e2  +  r3e3)  -  +  Yr3)]wcosx 

-  (r2c3  -  r3e2)wsinx)  +  +  (r32  +  rj2)W2 

-  V3W3  *  (ar3  -  yr^s  +  [PO^  +  r^) 

-  e2(yr3  +  ar^lwcosx  -  (x-j^  "  ri£3^wsinX^ 

.  k{-r3r1W1  -  r3r2W2  .  (Xj2  .  x22)Wj  »  (gXj  -  «x2)J 
+  (y(r1e1  +  r^)  -  £3(0^  +  8r2)]a>cosx 

-  (rxe2  -  r2e1)wsinx} 
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1.7 


ANGULAR  MOMENTUM 


The  angular  momentum  of  the  gun  with  respect  to  a  nonrotating 
2 

observer  at  point  P  is 


Vi  ♦  3i2w2  ♦  £1/3 


Consequently,  by  Eq.  (B-3)  in  Appendix  B,  the  angular  momentum  of  the  gun 
with  respect  to  point  0  is 


H  =  iljWj  ♦  jI2W2  .  kI3W3  *  MR„  X  l>0 


(1.22) 


Similarly,  the  angular  momentum  of  the  projectile  about  point  0  is 


h  =  ai^w^  +  b*-2w2  +  C^3W3  +  mRl  x  R1 


(1.23) 


The  vectors  a,  6,  c  are  resolved  into  i,  £  components  by  the  equations, 


A  A  a.  A 

a  =  ia1  +  +  ka7 


A  A 


b  =  ibj  +  jb2  +  kb3 


A  A  A 

c  =  iCj  +  jc2  +  kc3 


Consequently, 


h  =  i(?tjijWj  +  b  j  i  2w2  +  C1*3W3^  +  ^a2*lWl  +  b2*2W2  +  C2*3W3^ 


+  k(a3i!W1  +  b3i2W2  +  C3i3W3)  +  mRl  X  R1 


(1.24) 


Consequently,  in  view  of  Eqs.  (1.3)  and  (1.22), 


2 

JH.  L.  Langhaar  and  A.  P.  Bovesi,  Engineering  Mechanics  -Dynamics  (U) ,  McGraw- 
Hill }  New  York,  1959y  Art.  15-10. 
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H  +  h  =  iCIjWj  +  a1i]v1  +  bji2w2  +  Cli3w5) 


+  j (1 2^2  +  a2^1Wl  +  ^2^2W2  +  C2*3W3^ 


♦  k(I3Wj  +  ajijWj  *  b3i2v.2  *  c3i3w3)  ♦  r  x  r  (1.25) 


The  term  r  x  r  is  expressed  in  the  form  (-)i  +  (-)j  +  (-)k  by 
Eq.  (1.21).  Consequently,  Eq.  (1.25)  represents  H  +  h  in  the  form 


A  A  A 


H  +  h  =  iA  +  jB  +  kC 


(1.26) 


By  the  law  of  conservation  of  angular  momentum,  H  +  h  is  a  constant  vector. 
However,  in  general,  (A,B,C)  are  not  constants,  since  the  gun  moves  and 
accordingly  (i,j,k)  are  time -dependent  vectors.  The  constant  vectors  (i*,j' 

A 

k')  may  be  introduced  by  means  of  the  equations, 

A  A  A  /N 

i  *  x-1i'  +  m,  j  '  +  n..k! 


j  =  +  m2j 1  +  n2k' 


k  =  £ji 1  +  m^j '  +  n^k’ 


Thus,  H  +  h  is  expressed  as  follows: 


H  +  h  =  +  Jl2B  +  S.jC)  +  j  *  (m^A  +  m2B  +  m^C)  +  k'^/V 


+  n2B  +  n^C) 


Accordingly, 


J^A  +  S^B  +  Z5C  =  Kj 


m.A  +  nuB  +  m,C  =  K„ 


lijA  +  n,,B  +  n3C  = 


(1.27 


where  Kj,  K^,  are  constants.  However,  Eq.  (1.27)  introduces  a  complica¬ 
tion,  since  (fl-.m^n^)  are  unknown  functions  of  time. 

A  great  simplification  occurs  if  the  system  is  initially  at  rest, 
since  then  H  +  h  =  0  and  ^  =  K2  =  K3  =  °*  Consequently,  A  =  B  =  C  =  0. 
Attention  is  restricted  to  this  case. 

Provisionally,  the  following  notation  is  introduced: 


Mm  - 
M  +  m  r 


x 


r 


iX  +  jY  +  kZ 


(1.28) 


The  terms  X,  Y,  Z  are  given  by  Eq.  (1.21).  With  Eqs.  (1.25)  and  (1.28),  the 
law  of  conservation  of  angular  momentum,  H  +  h  =  0,  yields 


Vl  *  Vl“l  +  bl*2*2  +  V3"3  *  X  =  “ 


^2^2  +  a2^1Wl  +  ^2'*'2W2  +  C2^5W3  +  Y  =  0 


l?3  *  Vl“l  *  b3i2“2  *  V.*  *  Z  =  0 


(1.29) 


Now  Wj,  w2>  w^  are  given  by  Eq.  (1.6),  and  X,  Y,  Z  are  given  by  Eq.  (1.21). 
For  brevity,  the  following  notation  is  introduced: 


P  = 


Mm 


M  +  m 


(1.30) 


Thus,  the  following  equations  are  obtained  from  Eq.  (1.29): 

f1!  +  +  6j2i2  +  cx2i3  +  y(r22  +  r^2)^  +  [a^ij 


+  bxb2i2  +  Cj^ij  -  yr1r2]W2  +  [a^ij  +  b^^  +  CjC^ 


-  yi^r 


3]W2 


-a1i1(aa1  +  8a2  +  Ya3)u>  -  b^fabj  +  8b 2  +  yb^oo 

21 


-  c1i3(ae1  +  0c2  +  yc^w  -  u(yr2  -  $r^)§ 

-  iJa(r2e2  +  r3£3)wcosx  +  +  yr^wcosx 

+  P(r2e3  -  r3c2)wsinx 

Iala2i i  +  bib2i2  +  clc2i3  "  ^rir2^Wi  +  tX2  +  a22il  +  b22i2 
2  2  2 

+  c2  i3  +  y(r3  +  rj  ))W2  +  [a2a3i1  +  b2b3i2  +  C2C3i3  “  U12r3-lW3 
=  -a2i1(aa1  +  0a2  +  ya3)w  -  b^fcibj  +  0b2  +  yb^w 

-  ^ijCaCj  +  0c2  +  yc3)w  -  u(otr3  -  yr,)s  -  ii0(r3e3 

+  l'je^wcosx  +  Pe2(-Yr3  +  “rpwcosx  +  ‘  i^e^wsinx 

[a^ajij  +  b3b1i2  +  CjCjij  ■  Mr3r1]W1  +  [a^ij  +  b2b3i2 

+  C2C3i3  "  yr2r3^W2  +  ^3  +  a32il  +  +  C32i3 

2  2 

+  uOj  +  r2  )]W3  =  +  0a2  +  ya^w  -  b^i^ort^  +  0b2 

+  yb3)w  -  c3i3(ac1  +  0c,  +  yc3)a) 

-  3(0r1  -  ar2)s  --  nyCr^  +  r2e2)u)cosx 

+  pe3(ar1  +  0r2)wcosx  +  “  r^^wsinx 

(1.31) 

The  unknowns  in  Eq.  (1.31)  arc  W^,  W2>  W  .  Tlie  quantities  Ij, 

2,  Ij,  i^,  i2>  i3>  u,  a,  0,  y,  e^,  e 2>  e3  are  known  constants.  The 
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quantities  s,  x>  w  are  regarded  as  known  functions  of  t.  Perhaps  the 
simplest  way  to  determine  x  is  by  the  relationship. 


X  * 


(1.32) 


Since  ej,  e2,  are  known  constants,  Eq.  (1.1?)  determines  rJ(  r2,  r3  as 
functions  of  t.  By  Eq.  (1.10),  the  (i,j,k)  components  (a. ,b.,c.)  of 

/\  A  A  111 

vectors  (a,b,c)  are  known  functions  of  t.  With  a  computer  program,  Eq. 
(1.31)  can  be  solved  for  any  sequence  of  values  of  t  that  covers  the  period 
in  which  the  projectile  is  in  the  tube.  Thus,  the  functions  W^t),  W2(t), 
Wj(t)  can  be  tabulated  and  plotted. 

1.8  MOTION  OF  THE  GUN 

A  knowledge  of  the  functions  WJ(t),  W2(t),  W^(t)  does  not 
immediately  determine  the  motion  of  the  gun  in  reference  frame  F.  It  is 
necessary  to  determine  the  direction  cosines  (^i,m^,ni)  as  functions  of  t. 
When  these  functions  are  known,  the  motion  of  the  center  of  mass  of  the  gun 
is  determined  by  Eqs,  (1,4)  and  (1.17),  and  the  absolute  orientations  of  the 
principal  axes  of  the  gun  are  determined  as  functions  of  t. 

'Hie  problem  of  determining  the  functions  ft. ,  m.  ,  n.  is  purely  one 
of  kinematics  of  a  rigid  body.  The  results  are 


L  =  Wi  -  W  ft  ;  ft  =  W.ft  -  W  ft  ;  ft,  =  W^ft,  -  W  ft  ; 

1  o  Z  Z  o  Z  13  31  3  21  12* 

"*1  =  W3m2  ~  W2m3  ;  n*2  =  Wim3  '  W3mi  ;  ^3  =  W2ml  "  Wim->  ; 

**1  =  W3n2  ’  W2n3  ;  "2  =  Wln3  '  Vl  J  "3  =  W2nl  "  Wl”2  (1‘33) 


Since  the  functions  Wj,  W2,  W3  have  been  determined,  Eqs.  (1.33)  are 
differential  equations  th; 
initial  values  are  given. 


differential  equations  that  determine  the  functions  ft^,  nn,  when  the 


'N.  L.  Langhaar,  Energy  Methods  in  Applied  Mechanics  (V) ,  John  Wiley  and 
Sons,  Neu)  York,  W OP,,  Chap.  7 . 
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Instead  of  working  with  the  nine  unknown  functions  ,  m,  ,  n 

it  is  possible  to  work  with  three  liuler  angles  (0,i|>,t|O  by  means  of  the 

3 

following  relations  : 

-  sintjisin^  -  cos0cos<)>cosijj 

liij  =  -eos<j)siinjj  -  cos0sin<|>cosij> 

Uj  =  sin0cos(]> 

=  sin<f>cosi|>  +  cos0cos(|>sin«|> 

nu  =  -cos(f>cosiJ)  +  c.  os  Osinas  in 

n?  =  -sinOsinif) 

=  sinOeos(b 

A 

m„  -■  sinOsiiuJ) 

=  cosO 

The  differential  equations  for  0,  <f>,  <|>  ore'' 

•  • 

IVj  =  -OsiniJ)  +  (J)sin0cosi|> 

•  • 

W,,  =  -Ocosiji  -  (|>sinGs inij> 

•  • 

W„  =  'fcosO  +  tj) 


(1.34) 


( 1 . SS) 


2<\ 


Equations  (1.35)  are  easily  solved  algebraically  for  the  derivatives.  Thus, 


=  -W,  si  nil)  -  W.,cosi|) 
dt  l  l 

=  (W^cosij)  -  W^sini)))  cscO 


-  (W^cosi/<  -  W2sini|))cot0 


(1.36) 


According  to  an  existence  theorem  in  the  theory  of  ordinary  first 
4 

order  differential  equations,  there  is  a  unique  solution  0(t),  <J>(t. ) ,  'J'(t) 
of  Eq.  (1.36)  which  takes  given  initial  values  (9q,<|>0>i|j0)  ,  provided  that  9 
avoids  the  singular  values,  9  =  0  and  6  =  ir.  After  0(t),  <J>(t),  i|)(t)  are 
determined,  (£.,nu,ic)  are  determined  by  Eq.  (1.34). 

The  existence  theorem  also  applies  directly  to  F.q.  (1.33),  i.e., 
there  are  unique  functions  <L(t),  rn^(t),  iu(t)  that  satisfy  liq,  (1.33)  and 
that  take  given  initial  valuos  (It.  ^n  0) ,  Although  Eq.  (1.33)  involves 
nine  dependent,  variables,  it  has  the  advantage  that  the  equations  are 
linear. 

A  numerical  solution  of  the  differential  equations  appears 
fcasiblo.  It.  is  not  necessary  to  project  far  into  the  future,  since  the 
projectile  quickly  leaves  the  muzzle. 

1.9  CASE  OP  A  BALANCED  PROJECTILE 

If  tiie  projectile  is  perfectly  balanced  i  ^  =  i?  and  e,  =  c2  = 

=  0.  Also  Cj  =  a,  c2  =  0,  -  y.  Then  the  first  of  Eqs.  (1.31) 

becomes 

,T1  +  '|(al2  +  hl2)  *  0,2 ' 5  *  U(r22  +  l32)1Wl 
+  (i1(a1a2  +  bjb,)  +  aOi3  -  Mrlr2lW2 


*  MV3  *  l>ll>5>  +  "Yi3  -  '"'|''3'K3 


-I,. 


A’.  L.  TnaCy  Ovdinavij  DLj'j'eivnlLul  Equation a  (ID .  Doom'  Pubs.,  i\ln.i  Yovk ,  1944, 
A and  3.3/. 


2  2 

■ijWjaaj  +  +  yap3  +  abj  +  Bb^  +  ybpp 


2  2  2  • 
ai3(a  +8  +  Y  )w  -  p(yr2  -  Brps 


In  view  of  identities  among  the  direction  cosines,  this  reduces  to 
[Ix  +  ij  +  a2(i3  -  ip  +  y(r22  +  r32)]Wj  +  [<*B(i3  -  ip 


-  yr1r2]W2  +  [ay(i3  -  ip  -  yr^]!^  =  -  y(yr2  -  Br^s 


(1.37) 


Likewise  from  the  second  and  third  of  Eqs .  (1.31), 


[a3 C i 3  -  ij)  ~  yrJr2]W1  +  [I2  +  ij  +  B2(i3  -  ip 


+  ^r/  +  ri2)]W2  +  [BY(i3  '  V  -  Ur2r3]W3 


=  -BijW  -  U(ar3  -  Yrj)s 


[ay(i3  -  ip  -  ur^pWj  +  [6y(i3  -  ip  -  Ur3r2lW2 


+  lJ3  +  i\  +  Y2(i3  -  ij)  +  MCrj2  +  r22)]W3 


=  -YijW  -  h(Br1  -  ar2)s 


(1.37) 


(1.37) 


Equation  (1.37)  is  a  simplified  form  of  Eq.  (1.31)  that  applies  only  if  the 
projectile  is  perfectly  balance.  . 

l.io  Elementary  relations  between  forces  and  moments 

The  preceding  theory  deals  with  the  motion  of  a  free  rigid  gun. 

The  forces  and  moments  of  interaction  between  the  projectile  and  the  gun  also 
are  of  interest.  Although  gravity  has  been  neglected,  it  would  have 
only  a  small  effect  on  these  forces  and  moments. 


The  resultant  force  of  contact  that  the  projectile  exerts  on  the 
tube  is  designated  as  F.  The  force  of  contact  that  the  tube  exerts  on  the 
projectile  is  -F.  if  the  gun  is  regarded  as  a  free  body,  the  net  external 
force  on  it  is 

F  -  ApQV  Cl. 38) 

If  the  projectile  is  regarded  as  a  free  body,  the  net  external  force  on  it 
(neglecting  resistance  of  air  ahead  of  the  projectile)  is 


-F  +  ApjV  (1.39) 

The  detailed  forces  of  contact  that  the  projectile  exerts  on  the 
tube  are  designated  as  fj,  •  •  Hence, 

F  =  If,  (1.40) 

Let  force  f^  act  at  the  point  r  +  A^,  where  r  is  the  vector  from  the  center 
of  mass  of  the  gur.  to  the  center  of  mass  of  the  projectile.  Then  the 
moment  that  the  forces  f^  exert  about  the  center  of  mass  of  the  gun  is 

M  =  I(r  +  X.)  x  f.  =rxF+  IX.  x  f.  (a) 

v  l  i  ii 

The  forces  of  contact  that  the  tube  exerts  on  the  projectile  are  -fj,  -f ^ , 
-f^,  .  .  .  The  moment  of  these  forces  about  the  center  of  mass  of  the 
projectile  is 


M’ 


(b) 


By  Eqs .  (a)  and  (b) , 

M'  =  r  x  F  -  M 


(1.41) 
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The  moment  about  the  center  of  mass  of  the  gun  of  all  external 
forces  that  act  on  the  gun  is 


M  -  e  x  vp.A 


(1.42 


where  e  is  the  vector  from  the  center  of  mass  of  the  gun  to  the  initial 
location  of  the  geometric  center  of  the  projectile.  The  moment  about  the 
center  of  mass  of  the  projectile  of  all  the  forces  that  act  on  the  projectile 


M’  -  G  x  vp.A 


(1.43 


where  c  is  the  vector  from  the  geometric  center  of  the  projectile  to  the 
center  of  mass  of  the  projectile.  If  the  center  of  mass  of  the  projectile 
is  on  the  axis  of  the  tube,  e  =  0. 

1.11  MOMENT  OF  FORCES  ON  THE  GUN 

The  resultant  moment  M  exerted  on  the  gun  by  the  forces  of 
contact  with  the  projectile  is  resolved  into  components  along  the  principal 
axes  of  inertia  of  the  gun  through  its  center  of  mass,  i.e., 


M  =  iMx  +  jM2  +  kMj 


Also,  by  definition. 


A  a  A  /v 

v  -  ia  +  j3  +  ky 


_  A  A  A 

e  =  ie^  +  ie2  +  ke^ 


(1.44 


Consequently,  by  Eq.  (1.42),  the  net  moment  about  the  center  of  mass  of  the 
gun  is 


-  prtA(ye„  -  Be,)]  +  j  [M„  -  p„A(ae,  -  ye.)] 


The  components  M.  are  determined  by  Euler’s  dynamical  equations  for  a  rigid 

3  1 

body  ;  namely, 

dW. 

J1  ~dt“  ^2  "  VW2W3  +  P0A(Ye?  "  Se3-)  =  M1 

dw2 

l2  “dT  -  (I3  -  Wl  +  P0A(ae3  "  YeP  =  M2 


dW 

l3  It  ~  (I1  "  I2)W1W2  +  P0A(Bel  "  Cte2)  =  M3  (1-46) 

Since  (W^,W2,W3)  are  obtained  from  a  computer  program  for  the  solution  of 
Eq.  (1.31)  (or,  in  the  case  of  a  balanced  projectile,  from  Eq.  (1.37)),  the 
derivatives  dW./dt  can  be  obtained  by  numerical  differentiation.  The 
quantities  (a,0,y),  (ei»e2,e3^»  and  (Ij*^***)  are  known  constants. 
Consequently,  Eq.  (1.46)  determines  M,  provided  that  the  breech  pressure 
P0(t)  is  known. 

1.12  FORCE  OF  A  BALANCED  PROJECTILE  ON  THE  TUBE 

The  angular  acceleration  of  the  gun  is  W.  Since  W  =  iW  +  jW 

A  14 

+  kw3> 

*  A*  A*  A*  AAA 

W  =  iWL  +  jW2  +  kw3  +  iw:  +  jW2  +  kW3 

Also, 

i  =  W  x  i  =  jW3  -  kW2 

j  =  W  x  j  =  kWj  -  iW3 

k  =  W  x  k  =  iW2  - 

Consequently, 
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i.Wx  -  jW2  +  kW3  =  0 


Therefore, 


a  •  a  ,  A  • 

W  =  iWj  +  jW2  +  kW3 


(1.47) 


The  center  of  mass  of  the  gun  is  located  in  a  Galilean  reference 
frame  by  vector  R^,  which  is  given  by  Eq.  (1.2).  The  external  force  on  the 


gun  is  MRq .  Consequently,  by  Eqs .  (1.38), 


F  -  Apnv  =  -yr 


(1.48) 


where  y  =  Mm/ (M  +  m) . 

A  balanced  projectile  is  considered  in  this  article.  Consequently, 
by  Eq.  (1.13), 


r  =  e  +  vs 


(1.49) 


Therefore, 


rr  —  A  A  ( 

r  =  e  +  vs  +  vs 


However , 


e  =  W  x  e  and  v  =  W  x  v 


Therefore, 


r  =  W  x  r  +  vs 


(1.50) 


Differentiation  of  Eq.  (1.50)  yields 


r=Wxr+Wxr+vs+vs 


Hence, 


7  =  W  x  F  +  W  x  (W  x  r)  +  2W  x  vs  +  vs  (1 .51) 

By  expansion  of  the  vector  triple  product,  this  becomes 

r  =  W  x  F  +  W  W  •  F  -  Fw2  +  2W  x  vs  +  vs  (1 .52) 

Consequently,  by  Eq.  (1.48), 

F  a  Ap0v  -  u(W  x  r  +  W  W  •  r  -  rW2  +  2W  x  vs  +  vs)  (1.53) 

With  Eq.  (1.53),  the  net  force  on  the  gun  is  determined  by  Eq.  (1.38),  and 
the  net  force  on  the  projectile  is  determined  by  Eq.  (1.39).  In  view  of 
Eq.  (1.39),  the  axial  component  of  force  on  the  projectile  is 

-F  •  v  +  Apx  (1.54) 

For  numerical  computations,  Eq.  (1.53)  must  be  expressed  in 
scalar  form.  The  components  of  F  are  represented  by 

F  =  iFj  +  j  F2  +  kFj 


Also, 

A  A  /A  A 

v  =  ia  +  j8  +  ky 

and 


A  A  A 

r  =  ir1  +  jr2  +  kr3 
In  view  of  Eq.  (1.47), 

W  x  7  =  i(W2r3  -  W3r2)  +  jO*^  -  W^)  ♦  £(» ^  -  W^) 
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Accordingly,  Eq.  (1.53)  yields 


F1  =  a(Ap0  '  "  p^2r3  ‘  ^3r2)  “  WW1  (W2r2  +  W3r3) 

+  yr1(W22  +  W32)  -  2ys(yW2  -  BWj) 

F2  =  6(AP()  -  us)  -  ut^rj  -  Wjr3)  -  uW^W^  +  W3r3) 

4  Ur2(W32  4  Wx2)  -  2ys(aW3  -  yWj) 

F3  =  Y(Ap0  -  US)  -  uCW^r 2  -  W2rr)  -  uW3(W1r1  4  W2r2) 

4  ur3(W12  +  W22)  -  2ys($W1  -  aW2) 


(1.55 


The  solution  of  Eq.  (1.37)  provides  the  functions  W^,  W2>  Wj.  Consequently, 
the  force  F  is  determined  by  Eq.  (1.5S),  provided  that  the  center  of  mass 
of  the  projectile  lies  on  the  geometric  axis  of  the  tube.  Then  the  moment 
on  the  projectile  is  determined  by  Eq.  (1.41),  which,  in  expanded  form  is 


M  '  a  r  F  - 
2  3  1 


rlF3  “  M2 


M3*  =  rlF2  ”  r2Fl  M3  (1.56 

1.13  FORCE  OF  AN  UNBALANCED  PROJECTILE  ON  THE  TUBE 

If  the  projectile  is  dynamically  unbalanced,  r  and  r  are  given  by 
Eqs.  (1.17)  and  (1.19).  The  angle  x  through  which  the  projectile  has 
turned  relative  to  the  tube  is 


Differentiation  of  Eq.  (1.19)  yields 


^  ^  ^ 

r  ■  f  x  r  +  f  x  r  +  vs  ♦  vs  +  wv  x  c^cosy  +  wv  x  g^cosx 

•  • 

~  —  2A  —  •—  —  2— 

+  o>v  x  eQcosx  -  w  v  x  eQsinx  -  we0sn>x  *  we0sinx  -  w  eqcosx 

Now  r  can  be  eliminated  by  Eq.  (1.19).  Also, 

m  • 

A  _  A  _  _  _ 

v  =  W  x  v  and  =  W  x 

Consequently, 

r  =  Wx"r  +  Wx  [Wxr  +  vs+wOx  EgCosy  -  uiEgSinx] 

+  W  x  vs  +  Vs  +  wv  x  e"0cosx  +  u>(W  x  v)  x  eqcosx 
+  uv  x  (is  x  e0)cosx  -  W2v  x  e0sinx  -  uEgSinx 
-  (flis  x  ?Qsinx  -  w2?0cosx 

Since 

(W  X  V)  X  Eg  +  V  X  (W  X  Eg)  =  W  X  (V  X  Eg) 

this  reduces  to 

r  =  Wxr  +  Wx  [W  x  r  +  2vs  +  2wv  x  Egcosx  -  2w£0sinx] 

+  vs  +  u)(v  x  EqCosx  -  EQsinx)  -  w2(v  x  £0sinx  +  eqcosx) 

(1.58) 

If  efl  =  0,  Eq.  (1.58)  reduces  to  Eq.  (1.52). 
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For  brevity,  vectors  u  and  v  are  defined  by 
v  x  epcosx  -  ^sinx  *  u  =  iUj  *  ju2  +  £u3 
v  x  e^sinx  *  eQcosx  =  v  =  iv1  ♦  jv2  +  kv^ 

-  A  A  A 

Also,  “  ie^  +  je2  +  ke^.  Consequently, 

U1  =  ^e3  "  Ye2^co'  ”  EjSinx 


I 

i 


U2  =  (Yel  “  «e3)cosX  -  G2sin)( 
u3  =  (ae2  -  Be^cosx  -  e3sinx 

V1  =  *^e3  '  ^e2)sinx  +  eicosX 
v2  =  (yGj  -  ae3)sinx  +  e2cosx 

v3  =  (ae2  -  Bej)sinx  +  e3cosx 
By  Eqs.  (1.48)  and  (1.58), 

-  A  -  —  — .  • 2 — 

F  =  ApQv  -  virQ  -  2wyW  x  u  -  yuu  +  yto  v 

where  rQ  is  the  value  of  r  when  =  0.  Hence,  7^  is  given  by  Eq 
Set 


F-  *  Ap0V  -  HrQ 

Then  F1  is  the  expression  that  was  obtained  for  F  in  the  case  of 
projectile  (Eq.  (1.53)).  Set 


(1.59) 


(1.60) 

(1.61) 

(1.62) 

(1.52)  . 

(1.63) 

balanced 
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(1.64) 


F  =  ?•  ♦  f  ;  f  =  ifx  ♦  jf2  +  kf3 


Tlten  f  is  the  correction  to  F*  to  account  for  unbalance  of  the  projectile. 
In  scalar  form,  F»  is  given  by  Eq.  (1.55).  By  Eqs.  (1.62),  (1.63),  and 
(1.64), 


f  =  u[-2wW  x  u  -  u)u  +  to  v] 

Hence, 

fj  =  jj[-2w(W2u3  -  WjU2)  -  wu1  + 

f2  =  m[-2co(W3u1  -  Wxu3)  -  wu2  +  w2v2] 

f3  =  p[-2w(W1u2  -  W2Uj)  -  wuj  +  w2v3] 


(1.65) 


(1.66) 


The  quantities  (f^f^fj)  are  the  corrections  to  be  added  to 
(Fi,F2>F3),  respectively,  in  Eq.  (1.55)  to  account  for  dynamic  unbalance  of 
the  projectile.  It  is  to  be  recalled,  however,  that  W  is  affected  to  some 
extent  by  unbalance  of  the  projectile.  Consequently,  (Wj.l^.Wj)  are  to  be 
computed  by  Eq.  (1.31)  rather  than  by  Eq.  (1.37).  On  the  other  hand,  the 
quantities  (W^,W2,W3),  determined  by  Eq,  (1.37)  are  to  be  used  in  conjunction 
with  Eq.  (1.55).  The  moment  M'  on  the  projectile  is  given  by  Eq.  (1.56)  in 
either  case.  Equation  (1.46),  which  gives  the  moment  M  of  the  contact 
forces  that  the  projectile  exerts  on  the  gun,  is  valid  whether  or  not  the 
projectile  is  dynamically  balanced.  It  is  to  be  noted  that  Eqs.  (1.11)  and 
(1.59)  show  that  e  =  v. 
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SECTION  2 

FORCES  AND  MOMENTS  ON  A  RIGID  IMMOVABLE  GUN 
WITH  AN  UNBALANCED  PROJECTILE 


2 . 1  INTRODUCTION 

In  this  section,  the  gun  is  considered  to  be  rigid  and  immovable. 
The  motion  of  the  projectile  in  the  tube  is  prescribed.  The  net  force  and 
moment  acting  on  the  projectile  accordingly  are  determined  by  the  dynamical 
theory  of  a  single  rigid  body.  The  effect  of  gravity  is  neglected,  but  it 
would  merely  augment  the  force  F  by  the  term  mg,  where  is  the  vectorial 
acceleration  of  gravity.  Gravity  would  not  alter  the  moment  equations. 

2.2  KINEMATIC  RELATIONS 

Since  the  angular  velocity  of  the  gun  is  zero,  the  angular 
velocity  of  the  projectile  is 

w  =  vo)  (2.1) 


where,  as  before,  w  is  the  spin  of  the  projectile,  and  v  is  a  constant  unit 
vector  along  the  axis  of  the  gun  tube. 

Equations  (1.7),  (1.8),  (1.9),  (1.10),  and  (1.11)  again  apply. 
Since  the  origi »  0  of  coordinates  (x,y,z)  is  now  arbitrary,  it 
is  conveniently  taken  to  b°  the  initial  geometric  center  Q^'  of  the 
projectile  (Figure  1).  Then  Rq  =  e  =  0,  since  the  center  of  mass  of  the  gun 
is  irrelevant.  Also,  the  axes  (x,y,z)  may  be  oriented  in  the  directions 
aQ,  bQ,  cQ  of  the  initial  principal  'axes  of  the  projectile.  Then  aQ  =  i, 

A  A  A  A 

b^  =  j ,  and  =  k.  Consequently, 
ax  =  1  ;  31  =  0  ;  yl  =  0 


a2  =  0  ;  e2  =  1  ;  Y2  =  0 

a3  =  0  ’  ^3  =  0  ’  Y3  =  1 


(2.2) 
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In  view  of  liq.  (1.8)  and  Eq.  (2.1),  Wj  =  w  •  a  *  wv  •  a  =  wv  •  aQ  =  wv  •  i 
=  wet,  and  likewise  for  w2  and  w^.  Consequently, 

w1  =  toa  ;  w2  =  wB  ;  =  toy  (2.3) 

Hence , 

Wj  »  coa  ;  w2  =  uB  ;  w3  =  wy  (2.4) 

Since  v  and  are  constant  vectors  and  w  =  x,  by  Eqs .  (1.11)  and 

(1.59) 


e 


2- 
-(jo  e 


+  0)U 


(2.5) 


Since  Rg  =  e  =  0,  the  location  R^  of  the  center  of  mass  of  the 
projectile  is  (Figure  1) 


_  A 

r  =  Rj  *  vs  +  e 


(2.6) 


Consequently,  since  v  is  a  constant  vector,  the  acceleration  of  the  center  of 
mass  of  the  projectile  is 


r  =  vs  +  e 

Equations  (2.5)  and  (2.7)  yiold 

-  A..  2- 

r  =  vs  -  to  e  +  (jou 

With  Eq.  (1.10),  Eq .  (2.2)  yields 

a  a  2  2  a 

a  =  i  (a  +  (B~  ■>  y  )cosx]  +  j  [aB(l  -  cosx)  +  YsinxJ 

A 

♦  k[ya(l  -  cosy)  -  Bsiny] 


(2.7) 


(2.8) 


37 


b  =  i[od3(l  -  cosx)  -  ysinxJ  +  j  ifl  +  (Y“  *  ct')cosx] 


+  k  [3y(1  *  cosx)  +  asinxl 


c  =  i  [yot ( 1  -  cosx)  +  0si»x]  +  j  [3y ( 1  -  cosy)  -  «si»xl 


+  k  [Y2  +  (a2  +  B2)cosx3 


Equation  (1.11)  yields 


e  x  v  =  v)  x  (en  x  v)si»x  +  ert  x  vcosx 


Hence , 


(2.9) 


e  x  v  =  (en  -  vv  ♦  en)sinx  +  en  x  \>cosx 


(2.10 


2.3  FORCES  ON  THE  GUN  AND  ON  THE  PROJECTILE 

The  force  that  the  projectile  exerts  on  the  gun  by  direct  contact 
is  denoted  by  F,  as  before.  The  force  of  contact  that  the  gun  exerts  on  the 
projectile  is  -F.  The  net  force  on  the  projectile  is  accordingly 


-F  +  v(pjA  -  R) 


(2.11 


where  R  is  the  resistance  of  air  ahead  of  the  projectile.  In  view  of  F.q. 
(2.8),  Newton's  law  yields 


-F  +  v(p^A  -  R)  *  vms  -  mus  e  +  mum 


(2.12; 


Hence,  by  Cqs .  (1.11)  and  (1.59) 


F  =  v(PjA  -  R  -  ms)  +  mw  (v  x  "f^sinx  +  efgCosx) 


-  mw(v  x  e^cosx  -  ertsinx) 


(2.13 


The  axial  frictional  force  that  the  gun  exerts  on  the  projectile 

A  -  A  ,  A 

is  in  the  direction  -v.  Its  magnitude  is  F  •  v  =  1*^.  Since  \>  1  =  0, 

Eq.  (2.13)  yields 

If  =  P1A  ~  11  "  (2.14) 

By  definition, 

v  =  ia  +  +  icy  ,  3  iCj  +  je2  +  ke3  (2. IS) 

Consequently,  in  scalar  form,  Eq.  (2.13)  is 

Fj  =  a(p^A  -  11  -  ms)  +  mto2  [((k.j  -  yc2)  sinx  +  e^osyl 

-  mo  [(3g3  -  ye2)cosx  -  CjSinx] 

F2  =  3(p,A  -  H  -  ms)  +■  mw  [(ye^  -  ae3)sinx  +  e2cosx] 

-  mo[(yc1  -  oe.^)cosx  -  e2sinx] 

F3  =  y(p1A  -  K  -  ms)  •»  mo2((oie.2  -  Be^siny  +  GjCosx] 

-  mw[(oe2  -  Btjjcosx  -  e^sinxl  (2.16) 

The  vectors  and  v  x  cQ  arc  perpendicular  to  the  axis  of  the 
tube.  The  second  expression  in  Eq.  (2,15)  represents  centrifugal  forco.  It 
is  zero  if  the  eccentricity  is  zero. 

In  addition  to  the  force  F,  the  gun  experiences  the  breech  force 

A 

-vp^A  from  the  gases . 

2.4  MOMENT  ON  THE  PROJECTILE 

In  view  of  Eq.  (1.45),  the  moment  about  the  center  of  mass  of  the 

_  A 

projectile  of  all  the  forcos  that  act  on  the  pro  jec  tile  is  M*  -exv(p  A  -  .<) , 
(In  Eq,  (1.45),  R  was  neglected.)  The  components  of  this  moment  on  the 
principal  axes  of  the  inertia  of  the  projectile  are  obtained  by  taking  the 
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scalar  products  of  the  moment  with  a,  b,  and  c.  By  virtue  of  Eqs.  (2.3) 
and  (2.4),  Euler’s  equations  (Ref.  3)  yield 

M'  •  a  =  (PjA  -  R)e  x  v  •  a  +  ijCtw  -  (i2  -  i3)8yu)2 

M'  •  b  =  (pjA  -  R)e  x  v  •  b  +  i2Bw  -  (i3  -  ij)yotto2 

M'  •  c  =  (PjA  -  R)e  x  0  *  c  +  i3YW  -  (ij  -  i^aBui2  (2.17) 

If  U  is  any  vector,  aa  •  U  +  bb  •  U  +  cc  •  U  ■  U.  Consequently,  Eq.  (2.17) 
yields 

M'  -  (p^  -  R)e  x  V  +  ^(ijOia  +  i20b  +  ijYc)  -  u>2[(i2  -  i3)3Ya 

+  (i3  -  i^yab  +  (ij  -  i2)a$c]  (2.18) 

The  constant  vectors  v  ■  ict  +  j  6  +  ky  and  eQ  =  ie^  +  je2  +  ke7  are 
considered  to  be  known.  Also,  y(t)»  w  -  X»  Pi (t)  and  R(t)  are  regarded  as 

AAA 

known  functions.  Accordingly,  in  view  of  Eqs.  (1.11)  and  (2.9),  a,  b,  c, 

and  e  are  known  vector  functions  of  t.  Consequently,  Eq.  (2.18)  is  an 

explicit  vector  formula  for  the  moment  M’(t).  Equation  (2.10)  is  a 
representation  of  the  vector  product  e  x  v  that  occurs  in  Eq.  (2.18). 
Accordingly,  the  vector  M' ,  defined  by  Eq.  (2,18),  may  be  resolved  into  its 
i,  i,  k  components.  Setting  M'  =  iMj '  +  jM2'  +  kMj',  we  get 

Ml'  =  *-plA  "  R)  [EjSinx  -  a(ae1  +  0e2  +  ye3)sinx 

*  2  2 
+  (ye2  -  0e3)cosx]  +  [wijd  -  u'  (i2  -  i3)BY]  [« 

+  (32  +  y2)cosx]  +  [wi08  -  w2(i3  -  ip  va]  [aBversy 

A 

-  ysinx]  +  [wi3Y  -  w  (ij  -  i2)ag]  [yaversy  +  Bsiny] 
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m2*  =  (PjA  -  R)  [e2sinx  -  3(0^  +  3e2  +  ye3)sinx  +  (ae3  -ye^cosy] 
+  [wija  -  a)2(i2  -  i3)8y]  [aBversx  +  ysinx]  +  [d>i23 

-  to2(i3  -  ij)ya]  [3 2  +  (y2  +  a2)cosx]  +  [d>i3y 

2 

-  ii)  (ij  -  i2)aB]  [Byversx  -  otsinx] 


M3'  =  (pxA  -  R)  [e3sinx  -  yfaCj  +  3e2  +  ye3)sinx  +  (BEj  -  ae2)cosx] 
+  [arijOt  -  ai  (i2  -  i3)$y]  [yaversx  -  3sinx]  +  [wi23 
-  a)2 (i3  -  i1)ya]  [Byversx  •*  asinx]  +  [d)i3y  -  oj2 (ij  -  i2)aB] 
[y2  +  (a2  +  $2)cosx]  (2.19) 


where  versx  =  1  -  cosx. 

The  moment  M  of  the  contact  forces  that  the  projectile  exerts  on 
the  gun  about  the  origin  Q^'  is  given  by  Eq.  (1.41);  namely. 


M  =  r  x  F  -  M' 


(2.20) 


where  r  =  vs  +  e  . 

_  A 

If  the  projectile  is  balanced,  e  =  0  and  v  lies  along  a 

^  A  A 

principal  axis  of  inertia  of  the  projectile  -  say  v  =  k.  Then  Eq.  (2.19) 
yields  Mj'  =  M2’  =  0  and  Mj'  -  ijU).  In  this  case,  M3'  is  the  rifling 
torque.  For  an  unbalanced  projectile,  the  rifling  torque  may  be  defined  as 
M'  •  v,  where  M'  is  represented  by  Eq.  (2.18).  The  term  (pjA  -  R)£  x  0 
cancels  out  of  the  scalar  product  M'  •  v. 
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SECTION  3 

FORCES  AND  MOMENTS  ON  A  FREELY-RECOILING 
RIGID  GUN  THAT  IS  CONSTRAINED  AGAINST  ROTATION 


3.1  INTRODUCTION 

A  rigid  gun  translates  freely  along  an  inclined  guide,  represented 
by  an  inclined  plane  in  Figure  2.  It  is  constrained  against  rotation.  The 
angle  0  of  the  inclined  plane  may  differ  from  the  angle  of  elevation  0  +  <J> 
of  the  barrel.  Air  resistance  to  motion  of  the  projectile  and  effects  of 
gravity  are  included  in  the  analysis.  It.  is  questionable  whether  gravity  is 
meaningful  in  this  problem,  since  it  would  cause  the  system  to  slide  down 
the  inclined  plane  with  increasing  speed.  It  is  easily  eliminated,  however, 
by  setting  g  =  0. 

In  addition  to  notations  introduced  previously,  a  few  new 
notations  are  added. 

A 

U  is  a  unit  vector  along  the  axis  of  the  recoil  guide  (Figure  2) . 
£,n  are  rectangular  coordinates  fixed  in  a  Galilean  reference 
frame  that  contains  the  recoil  guide  (Figure  2) . 

u  is  the  recoil  displacement  along  the  guide  (Figure  2) . 

is  the  vector  from  the  origin  of  the  coordinates  (£,n)  to  the 
center  of  mass  of  the  projectile  (Figure  2) . 

3.2  MOTION  OF  THE  SYSTEM 

The  velocity  of  the  projectile  relative  to  the  gun  is  s.  The 
component  of  this  velocity  along  the  guide  is  scos<(»,  as  is  seen  by  Figure  2. 
Consequently,  the  momentum  of  the  system  ijr»  the  direction  of  the  guide  is 

•  •  •  t  • 

m(scos<}>  -  u)  -  Mu  =  mscos<f>  -  (M  +  m)u 

The  component  of  external  force  on  the  system  in  the  direction  of 
the  guide  is 

-  (M  +  m)gsin8  +  (Ap^  -  Ap^  -  R)cos<}> 

where  R  is  the  resisting  force  of  air  in  the  tube  ahead  of  the  projectile. 
Since  the  external  force  is  equal  to  the  rate  of  change  of  momentum, 
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Figure  2.  Gun  Schematic  and  Notations 


mscos<j>  -  (M  +  m)u  =  -(M  +  m)gsin0  +  (Ap^  -  Ap^  -  R)cos<f> 


or 


u  =  gsinQ  + 


m 


M  +  m 


(Apj  -  Ap  -  R) 

'icos* - sm — cos* 


(3.1) 


Since  the  gun  is  constrained  against  rotation,  the  absolute  angu¬ 
lar  velocity  of  the  projectile  is 


_  A 

w  =  v« 


(3.2) 


Since  to,  s,  and  R  are  regarded  as  known  functions  of  t,  Eqs.  (3.1)  and  (3.2) 
determine  the  motion  of  the  gun  and  the  projectile. 

3.3  FORCES  IN  THE  SYSTEM 

The  resultant  force  of  contact  that  the  projectile  exerts  on  the 
tube  is  F.  The  resultant  force  of  contact  that  the  tube  exerts  on  the 
projectile  is  -F.  If  the  gun  is  regarded  as  a  free  body,  the  net  force  on 
it  is 


F  -  Ap  v  +  Mg 


(3.3) 


where  p^  is  the  gas  pressure  at  the  breech.  The  net  force  on  the  projectile 


is 


-F  +  mg  +  (Ap^  -  R) v 


(3.4) 


where  Pj  is  the  gas  pressure  on  the  base  of  the  projectile.  Consequently, 
by  Newton's  law 


•F  +  mg  +  (Ap^  -  R)v  =  mRj 


(3.5) 


By  Figure  2, 


=  e  +  Vs  +  e 


(3.6) 


4 


Also, 


/V 


where  e^  is  the  initial  value  of  e.  Hence, 

_  _  A  A  - - 

R^  =  eQ  -  yu  +  vs  +  e 


(3.7) 


(3.8) 


Therefore, 


Rj  =  -mu  +  vs  +  e  (3.9) 

Equations  (1.11)  and  (2.5)  are  again  applicable.  Equations  (3.1),  (3.5), 
(3.9),  (1.11),  (1.59),  and  (2.5)  yield 

A  ^  2  A 

F  =  mg  +  (APl  -  R)v  +  ymgsin9  +  scos<f>  -  ^  H\(APi  “  Ap0 

2  A  —  _  *  A  — 

-  R)cos<{>  -  mvs  +  mu>  (v  x  eQsinx  +  e0cosx)  -  mw(v  x  eqcosx 

-  eQsinx)  (3.10) 

-  A 

The  magnitude  of  the  axial  friction  force  on  the  projectile  is  F  •  v  =  Ff. 
Hence, 

.  2,  mAp 

Ff  =  -mgcos0sin<j>  +  (APl  -  R  -  ms)  — ^  — -  cos  <j>  (3.11) 

If  M  °°,  F^  reduces  to  Eq.  (2.14),  aside  from  the  g-term  which  was  neglected 
in  the  derivation  of  Eq.  (2.14). 

3.4  MOMENT  ON  THE  PROJECTILE 

Since  the  motion  of  the  projectile  relative  to  the  gun  is 
presumed  to  be  prescribed,  the  recoil  of  the  gun  merely  superimposes  a 
translation  on  the  absolute  motion  of  the  projectile.  The  recoil  has  no 
effect  on  the  angular  velocity  of  the  projectile.  Therefore,  in  view  of 
Euler's  equations  (Ref.  3) ,  it  has  no  effect  on  the  resultant  moment  about 
the  center  of  mass  of  the  projectile.  Likewise,  gravity  has  no  effect  on 
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this  moment.  The  theory  of  moments  on  a  projectile  in  a  motionless  rigid  gun 
consequently  is  directly  applicable  to  the  recoiling  gun.  The  moment  M* 
about  the  center  of  mass  of  the  projectile  is  again  given  by  Bq.  (2.18). 

The  components  of  M*  are  again  given  by  Eq.  (2.19). 
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•  H  * 


SECTION  4 

RECOILING  RIGID  GUN  WITH  OFFSET  BREECH 
AND  FIXED  TRUNNION 


4.1  INTRODUCTION 

If  the  center  of  mass  of  the  breech  of  a  gun  lies  below  the  axis 
of  the  tube,  the  recoil  causes  the  muzzle  to  jerk  upward  when  the  gun  is 
fired.  Because  of  its  spin,  the  projectile  then  exerts  a  gyroscopic 
couple  that  tends  to  turn  the  tube  sideways.  It  is  assumed  in  this  section 
that  a  constraint  is  provided  which  prevents  rotation  of  the  gun  about  a 
vertical  axis.  Then  the  spin  of  the  projectile  has  no  effect  on  the  motion 
of  the  gun.  Also,  because  of  this  constraint,  offsetting  of  the  breech 
block  to  the  right  or  the  left  has  no  kinematic  effect.  When  the  gun  is 
fired,  each  particle  describes  a  curve  that  lies  in  a  plane  perpendicular  to 
the  axis  of  the  trunnion.  Because  of  the  constraint  provided  by  the 
trunnion,  the  gun  has  only  two  degrees  of  freedom.  The  projectile  adds 
another  degree  of  freedom  to  the  system. 

Figure  3  is  a  schematic  side  view  of  the  gun.  The  trunnion  is 
fixed  in  a  Galilean  reference  frame;  e.g.,  the  earth.  The  recoil  mechanism 
is  represented  schematically  as  a  spring  in  a  slot.  One  end  of  the  spring 
is  attached  to  the  breech  block,  and  the  other  end  to  the  trunnion.  The 
slot  slides  freely  over  the  fixed  trunnion.  The  spring  need  not  be  Hookean. 
Rather,  the  force  F  exerted  by  the  recoil  mechanism  is  regarded  as  an 
unspecified  function  of  u  and  u,  where  u  is  the  displacement  of  the  breech 
along  the  axis  of  the  slot  (Figure  3) .  Accordingly,  the  recoil  mechanism 
may  contain  nonlinear  springs  and  nonlinear  dashpots.  Also,  Coulomb 
friction  is  admitted.  For  generality,  the  line  of  action  of  the  recoil 
mechanism  is  not  taken  parallel  to  the  axis  of  the  tube. 

The  projectile  is  considered  to  be  a  body  of  revolution  with 
its  center  of  mass  on  its  axis  of  symmetry,  and  with  one  principal  axis  of 
inertia  coinciding  with  the  axis  of  symmetry  i.e.,  the  projectile  is 
perfectly  balanced. 

4.2  NOTATIONS 

Some  deviations  from  previous  notations  are  necessary.  Also,  a 
few  notations  are  added. 
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x,y  are  rectangular  coordinates  with  the  y-axis  vertical,  and  the 
origin  at  the  trunnion.  They  are  fixed  in  a  Galilean  reference  frame 
(figure  3) . 

ot  is  the  angle  between  the  axis  of  the  tube  and  the  axis  of  the 
recoil  mechanism.  It  is  a  constant  (Figure  3) . 

u  is  the  displacement  of  the  gun  along  the  axis  of  the  recoil 
mechanism  (Figure  3) . 

£,n  are  rectangular  coordinates  scribed  on  the  breech.  The  £-axis 
is  the  axis  of  the  recoil  mechanism. 

9  is  the  angle  between  the  £  and  x  axes.  Generalized  coordinates 

are  0,  u. 

CQ,n0  are  the  n  coordinates  of  the  center  of  mass  of  the 
projectile  before  firing. 

S,n  are  the  n  coordinates  of  the  center  of  mass  of  the  gun. 

x,y  are  the  x,  y  coordinates  of  the  center  of  mass  of  the  gun. 

Xq.Yq  are  the  x,  y  coordinates  of  the  center  of  mass  of  the 
projectile  before  firing. 

are  the  n  coordinates  of  the  center  of  mass  of  the 
projectile  at  time  t. 

xl’^l  are  t^e  x’  y  C00rclinates  °f  the  center  of  mass  of  the 
projectile  at  time  t. 

Fq  is  the  force  of  gas  pressure  on  the  base  of  the  tube  at  time  t. 

Fj  is  the  net  force  driving  the  projectile  at  time  t  (force  of  gas 

pressure  minus  the  friction  of  the  barrel) .  F1  includes  the  resistance  of 
air  ahead  of  the  projectile. 

F(u,u)  is  the  force  that  the  i*ecoil  mechanism  exerts  on  the 

breech. 

s  is  the  distance  that  the  projectile  has  traveled  relative  to 
the  gun  at  time  t. 

M  is  the  mass  of  the  gun. 

in  is  the  mass  of  the  projectile. 

I  is  the  moment  of  inertia  of  the  gun  about  a  transverse  axis 
through  the  center  of  mass  of  the  gun. 

i  is  the  moment  of  inertia  of  the  projectile  about  a  transverse 
axis  through  the  center  of  mass  of  the  projectile. 

T  is  the  kinetic  energy  of  the  system. 
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g  is  the  acceleration  of  gravity. 

are  components  of  generalized  force.  <5W  =  Q,6u 
+Q260+Q36s.  1 

4.3  COORDINATE  TRANSFORMATION 

By  Figure  4, 

x  =  (^  -  u) cosQ  -  nsinQ 

y  =  (5  -  u)sin0  +  rjcos0 


Consequently, 

x  =  (C  -  u)cos0  -  rjsinO 
y  -  (£  -  u)sinO  +  ncor>0 
x0  ~  -  u)cos0  -  nQsin0 

y0  =  (?0  -  u)sin6  +  nQcos0 

Xj  =  -  u)cos0  -  n^sin0 

yl  =  ^1  "  u)sinO  +  rijCosO 


Also,  by  Figure  3, 


^1  ”  ^“0  +  scosa 
Oj  =  Hq  +  ssina 


Therefore, 


x.  =  (Cft  -  u)cosO  -  rusinO  +  scos(0  +  a) 


^1  a  "’C  *  u)sm0  n^cosB  +  ssin(0  +  a) 


By  liq.  (4.2), 


•  • 

x2  +  y2  =  (f  -  u)202  +  (u  +  m0) 2 


By  Eq.  (4.4), 


Xj2  +  y j2  =  [O^  -  u)  +  nx]2  +  [^  -  u  -  rijO]2 


By  Eq.  (4.5), 


scosa  ; 


-  ssina 


KINETIC  ENERGY 

The  kinetic  energy  of  the  gun  is 


•  ♦ 

I  •  2  1  — ?  —9 

Tg  =  |  io' ♦  £  M(x“  +  y"> 


The  kinetic  energy  of  the  projectile  is 


1  .i.2  1  ,*  2  ♦  2. 

fp  =  2  10  +  2  ra(xl  +  yl  ) 


Consequently,  in  view  of  Eqs .  (4.5),  (4.7),  (4.8),  and  (4.9),  the  kinetic 
energy  of  the  system  is 


T  -  -kl  +  i)62  +  \  M[('C  -  u)202  +  (u  +  n6)2] 


+  jDi[0(^  +  scosa  -  u)  +  ssina]2 


\  •  ♦  •  2 

+  m  [scosa  -  u  -  0(nA  +  ssina)] 


4.5 


GENERALIZED  FORCE 

It;  the  coordinates  (u,G,s)  receive  virtual  increments  (6u,60,6s), 
the  work  of  all  the  forces  that  act  on  the  system  is  a  linear  form  in  those 
increments;  i.e., 

6W  =  Q^u  +  Q2<S0  +  QjSs  (4.11) 

Ql ,  Q2,  Q3  are  called  "components  of  generalized  force." 

The  virtual  work  of  gravity  is 

6wgr  =  'Mg6y  "  mg5yi 


By  Eq.  (4.2), 


6y  -  -(sin0)6u  +  f(‘C  -  u)cos0]S0  -  (nsinO)60 
By  Eq.  (4.6), 

6y^  =  KSq  -  u)cos0  -  nosin0  +  scos(0  +  a)]60  -  6usin0 
+  sin(0  +  a)6s 


Consequently, 

6Wr.  =  (M  +  m)gs.inO  6u  -  Mg  [(C  -  u)cosO  -  rfsinG]<S0 
gr 

-  mg((^Q  -  u)cos0  -  riQsin0]60  -  mgscos(0  +  oc)60 

-  mgsin(0  +  a) 6s  (4.12) 

The  virtual  work  of  gas  pressure  on  the  breech  (see  Figure  3)  is 
$Wbr  =  -Fqcos(0  +  ot)6x0  -  FQsin(e  +  «)6yQ 
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Consequently,  by  Eq.  (4.3), 


<SWbr  =  FgCOsaSu  -  Fq  [(Cq  -  u)sina  -  nQcosa]6e  (4.13) 

The  virtual  work  that  the  recoil  mechanism  performs  on  the  breech  is 

SW  =  -FSu  (4.14) 

The  component  of  absolute  virtual  displacement  of  the  projectile 
along  the  axis  of  the  tube  is 

6xjCos(9*a)  +  dy^sinfS  a) 

With  Eq.  (4.6),  this  reduces  to 

6s  -  6ucosa  +  (£q  -  u)sina69  -  n0cosa59 

Consequently,  the  virtual  work  performed  on  the  projectile  is 

6Wpr  =  F^s  -  F16ucosa  +  Fj(S0  -  u)sina69  -  F^coso^  (4.15) 

Equations  (4.11)  to  (4.15)  yield 

Qj  =  (M  +  m)gsin9  +  (FQ  -  F^cosa  -  F(u,u) 

Q2  =  -Mg[(C  -  u)cos0  -  nsin9]  -  mg[(£0  -  u)cos9  -  nQsin0] 

-  mgscos(9  +  a)  -  (FQ  -  Fj)  [(SQ  -  u)sina  -  n0cosa] 

Q  =  -mgsin(9  +  a)  +  F^  (4.16) 

3 
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The  component  may  be  identified  as  the  counterclockwise  moment  about  the 
trunnion  of  all  external  forces  acting  on  the  system.  The  forces  FQ  and 
are  regarded  as  external  forces . 

4 . 6  LAGRANGE ' S  EQUATIONS 

The  Lagrange  equations  are 


±c5j  _  IX  -  q 

dt  l9uJ  3u  ^1 


d  3T 

dtW 


dt  ^  3s 


(4.17) 


The  term  3T/30  is  missing  because  T  does  not  depend  on  0  (see  Eq.  (4.10)). 

The  derivative  3T/30  may  be  identified  as  the  angular  momentum  of  the  system 
about  the  trunnion.  Consequently,  the  second  Lagrange  equation  expresses  the 
fact  that  the  moment  of  all  the  forces  about  the  trunnion  equals  the  rate  of 
increase  of  the  angular  momentum  of  the  system  about  the  trunnion.  It 
might  be  argued  that  FQ  and  Fj  are  internal  forces,  and,  in  the  absence  of 
gravity,  the  ang:  ar  momentum  of  the  system  about  the  trunnion  is  constant. 
However,  this  is  true  only  if  the  gases  in  the  tube  are  included  in  the 
system.  Consequently,  it  is  best  to  regard  FQ  and  F1  as  external  forces. 
Force  F^  may  include  the  resistance  of  air  in  the  tube  ahead  of  the 
projectile,  which  is  clearly  an  external  force. 

With  Eq.  (4.16),  the  first  Lagrange  equation  yields, 

(M  +  m)u  +  (Mp  +  mrig  +  mssina)0  -  mscosa  -  (M  +  m)u0 

_  «  2  •  •  • 

+  (M£  +  m£g  +  mscosa) 0  +  2ms0sina  +  F(u,u) 


=  (M  +  m)gsin0  +  (Fg  -  F^)cosa 


5s 


(4.18) 


By  Eq.  (4.10), 


ft  T  __  •  ^  2  o 

—  =  [Mn  +  m(nn  +  ssina)  ]u  +  [I  +  i  +  M(£  -  u)  +  Mri 

30  u 

2  2  • 

+  m(£g  -  u  +  scosa)  +  m(Hg  +  ssina)  ]0 


+  m[(C0  -  u)sina  -  HgCosaJs 


(4.19) 


Accordingly,  by  Eqs.  (4.16)  and  (4.17),  the  second  Lagrange  equation  is 


S'. 

♦•v,* 


■t 


{  [Mn  +  m(nQ  +  ssina)]u  +  [I  +  i  +  M('C  -  u)2  +  Mn2 


+  m(£g  -  u  +  scosa)2  +  m(nQ  +  ssina) 2]0  +  m  [(£g  -  u)sina 


-  ngCosa]s)=  -  Mg[(£  -  u)cos0  -  nsin0]  -  mg[(£g  -  u)cos0 


-  ngSinO]  -  mgscos(0  +  a)  -  (FQ  - 


[(CQ  -  u)sina  -  nQcosa] 


(4.20) 


0 


The  Lagrange  equation  corresponding  to  s  is 


[(£g  -  u)sina  -  nocosa]0  +  s  -  iicosa  -  8  (£Q  -  u)cosa 


*2  •  •  *2  1 
-0s-  2u0sina  -  0  nQsina  =  -gsin(0  +  °0  +  ~jjj" 


(4.21) 


Equations  (4.18)  to  (4.21)  are  simplified  considerably  if  a  =  0. 

The  displacement  s  of  the  projectile  and  the  base  pressure  force 
Fg  may  be  regarded  as  known  functions  of  t.  Equation  (4.21)  may  be  used  to 
eliminate  from  Eqs.  (4.18)  and  (4.20).  After  is  eliminated,  Eqs. 
(4.18)  and  (4.20)  are  two  nonlinear  coupled  second-order  ordinary 
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differential  equations  that  determine  the  functions  u(t)  and  9(t),  if 
initial  values  uQ,  uQ,  8Q,  dQ  are  given.  After  u(t)  and  0(t)  are  determined 
F^(t)  can  be  calculated  by  Eq.  (4.21). 

Rotary  friction  and  rotary  spring  resistance  in  the  trunnion  have 
been  disregarded,  but  their  inclusion  in  the  equations  is  simple.  The 
right  side  of  Eq.  (4.20)  is  merely  augmented  by  a  term  -cj> (0,0)  ,  which 
represents  the  resisting  moment  of  the  trunnion. 

•  • 

If  the  system  starts  from  rest,  the  quantities  u,  u,  9  are  small 
while  the  projectile  remains  in  the  barrel.  Consequently,  it  is  reasonable 
to  linearize  the  differential  equations  in  these  variables.  At  least,  this 
approximation  provides  a  start  for  an  iterative  solution  of  the  nonlinear 
equations. 


SECTION  5 
CONCLUSIONS 


The  equations  are  rather  complicated,  but  they  appear  to  be 
amendable  to  numerical  treatment  with  a  digital  computer.  No  numerical 
results  are  included  in  this  report. 

Section  1  treats  a  gun  that  is  unsupported.  In  an  actual  gun, 
the  effect  of  the  recoil  mechanism  may  be  negligible  during  the  few 
milliseconds  that  the  moving  projectile  remains  in  the  barrel.  The  trunnion 
provides  a  constraint,  but  it  may  be  temporarily  ineffective  if  there  are 
appreciable  clearances  in  the  bearings.  Consequently,  the  unsupported  gun 
might  not  deviate  unduly  from  reality  in  some  cases. 

Equations  (1.31)  are  the  key  equations  in  Section  1.  For  a 
balanced  projectile,  they  reduce  to  Eq.  (1.37).  The  unknowns  in  Eq.  (1.31) 
or  (1.37)  are  the  angular  velocity  components  (W^^jW^).  Equations  (1.31) 
or  (1.37)  are  linear  algebraic  equations  in  these  variables.  Consequently, 
they  are  immediately  solvable.  After  (W. ,W2,W3)  are  calculated,  the 
direction  cosines  (JL,nu,n.)  of  the  principal  axes  of  inertia  of  the  gun 
are  determined  as  functions  of  time  t  by  solving  the  nine  linear  first- 
order  differential  equations  (1 .33)  or,  alternatively,  by  solving  the  three 
nonlinear  first-order  differential  equations  (1.36).  Apparently,  these 
equations  must  be  solved  numerically.  After  (JL,m.,m)  are  determined, 
the  motion  of  the  center  of  mass  is  determined  by  Eq.  (1.4).  Thus,  the 
motion  of  the  system  is  determined  completely,  since  the  motion  of  the 
projectile  relative  to  the  gun  is  presumed  to  be  known.  Since  the  angular 
velocity  components  (W^.W^.W^)  have  been  calculated,  the  components  (M^, 
N^.Mj)  of  the  moment  M  (exerted  about  the  center  of  mass  of  the  gun  by  the 
contact  forces  of  the  projectile)  are  determined  by  Euler's  dynamical 
equations  for  a  rigid  body  (Eq.  (1.46)).  The  force  F  that  a  balanced 
projectile  exerts  on  the  tube  is  determined  by  Eq.  (1.55).  For  an 
unbalanced  projectile,  this  force  must  be  augmented  by  the  corrective  terms 
in  Eq.  (1.66).  The  rifling  torque  and  the  axial  friction  force  on  the 
projectile  can  be  calculated  directly  after  the  vectors  M  and  F  are 
determined. 
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Section  2  treats  a  gun  that  is  immovable.  The  force  F  that  the 
projectile  exerts  on  the  gun  is  given,  in  this  case,  by  Eq.  (2.13),  or,  in 
scalar  form,  by  Eq.  (2.16).  The  moment  M*  about  the  center  of  mass  of  the 
projectile  of  the  contact  forces  imposed  by  the  tube  is  given  by  Eq.  (2,18), 
or,  in  scalar  form,  by  Eq.  (2.19).  These  equations  are  explicit  algebraic 
formulas  for  F  and  M' .  Other  pertinent  forces  and  moments  are  immediately 
determinate  when  F  and  M'  are  known. 

Section  3  treats  a  gun  that  translates  freely  along  a  guide,  but 
it  is  constrained  against  rotation.  In  this  case,  the  moment  M'  is  the 
same  as  for  the  fixed  gun,  treated  in  Section  2.  The  recoil  displacement  u 
is  determined  by  integrating  Eq.  (3.1).  The  force  of  contact  F  that  the 
projectile  exerts  on  the  tube  is  given  by  Eq.  (3.10).  The  magnitude  of  the 
axial  frictional  force  on  the  projectile  is  given  by  Eq.  (3.11). 

Section  4  treats  a  gun  with  a  trunnion  and  a  general  type  of 
recoil  device  (see  Figure  3) .  The  system  has  two  degrees  of  freedom, 
corresponding  to  the  recoil  displacement  u  and  the  angular  displacement  0 
of  the  gun.  A  third  coordinate  s  is  introduced.  It  represents  the  axial 
displacement  of  the  projectile  with  respect  to  the  gun,  but,  since  this  is 
presumed  to  be  given,  the  Lagrange  equation  corresponding  to  s  merely 
determines  the  net  driving  force  Fj  on  the  projectile.  After  Fj  is 
eliminated,  Eqs.  (4.18)  and  (4.20)  are  two  nonlinear  second-order  ordinary 
differential  equations  that  determine  the  functions  u(t)  and  0(t),  if 
initial  values  uQ,  u^,  8Q,  9^  are  given.  A  numerical  program  to  carry  out 
this  solution  is  needed. 

Within  the  frameworks  of  the  mathematical  models  that  are  used, 
the  analyses  are  exact.  However,  the  problem  of  gas  dynamics  in  the  tube  is 
not  rigorously  separable  from  the  problem  of  dynamical  response  of  the  gun. 
The  effects  of  momentum  and  kinetic  energy  of  the  charge  requires  further 
study. 

Finally,  it  may  be  advisable  to  issue  as  separate  complete  working 
reports  the  results  of  Sections  1,  2,  3,  and  4. 
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APPENDIX  A 

DISPLACEMENT  VECTOR  FIELD  OF  A  RIGID  BODY 

This  appendix  presents  a  derivation  of  the  displacement  vector 
field  of  a  rigid  body  that  undergoes  an  angular  displacement  x  about  an  axis 
A  with  direction  0  relative  to  given  rectangular  coordinates  (S.n.C).  The 
axis  A  is  conveniently  taken  to  pass  through  the  origin  0  (Figure  A-l). 

A  particle  P  of  the  body  describes  a  circular  arc  C  of  radius  a, 
whose  plane  is  perpendicular  to  axis  A,  and  whose  center  M  is  on  axis  A. 

The  displacement  vector  q  of  particle  P  may  be  resolved  into  components  q^ 

and  q9,  such  that  q1  is  tangent  to  circle  C  at  the  initial  point  P,  and  q 

is  parallel  to  the  radius  MP.  The  radius  vector  OP  is  denoted  by  p. 

By  Figure  A-l,  it  can  be  seen  that  |v  x  p|  =  a,  where  the 
v  x  p|  denotes  the  magnitude  of  the  vector  v  x  p.  Also,  v  x  p 
has  the  direction  of  q^ .  Consequently, 

qx  =  v  x  psinx 

The  magnitude  of  q9  is 
q2  =  a(l  -  cosx) 

—  a  a  —  .A  A  — -  * 

The  direction  of  q2  is  that  of  the  vector  v  x  (v  x  p) .  Also,  |v  x  (v  x  p) | 

a.  Therefore, 

q2  =  V  X  (V  X  P) (1  -  cosx) 

Since  the  displacement  vector  of  particle  P  is  q  =  qj  +  q2> 


q  =  v  x  psinx  +  v  x  (v  x  p)(l  -  cosx)  (A-l) 

The  vector  triple  product  in  Eq.  (A-l)  may  be  expanded  by  means 
of  the  identity 

Ax(BxC)  =  B  A  •  C  -  C  A  •  B 
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'Hie  particle  P  is  displaced  from  position  p  to  the  position  R  =  p  +  q. 
Consequently, 

R  =  v  x  psiny  +  pcosx  +  \>G  •  p(l  -  cosx)  (A-2) 
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APPENDIX  B 
MOMENTUM  PRINCIPLES 


An  arbitrary  mechanical  system  is  referred  to  a  Galilean  reference 
frame.  The  momentum  of  the  system  is 

G  =  [vdm 


where  v  is  the  velocity  of  the  mass  particle  dm,  and  the  integral  extends 
throughout  the  system  (Figure  B-l) .  Also,  v  =  dr/dt,  where  r  is  the  radius 
vector  from  the  origin  0  to  particle  dm  (Figure  B-l).  Consequently, 


Furthermore , 

jrdm  =  mrQ 

where  m  is  the  mass  of  the  entire  system,  and  rQ  is  the  radius  vector  from 
point  0  to  the  center  of  mass  of  the  system.  Therefore, 

d?0 

G  =  m  -jjr  =  mvQ  (B-l) 

where  v.  is  the  velocity  of  the  center  of  mass  of  the  system. 

U  2 
The  angular  momentum  of  the  system  about  point  0  is 

Hq  =  jr  x  vdm 

Likewise,  the  angular  momentum  of  the  system  about  another  point  P  is 
Hp  =  jp  x  vdm 

_ 

H.  L.  Langhaar  and  4.  P.  Bore  si.  Engineering  Mechanics-Uynamies  ( V) ,  MoGraio- 
llill ,  New  York,  1959,  Art .  15-10. 
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where  p  is  the  vector  from  point  P  to  the  mass  particle  dm.  By  Figure  B-l, 
r  =  D  4-  p  where  0  is  the  vector  OP.  Consequently, 


■v  ■  j 


(r  -  D)  x  vdm  -•  r  x  vdm  -  I)x  7dm 


Therefore, 


Hp  =  Hq  -  D  x  G 


Equations  (B-l)  and  (B-2)  yield 


Hp  =  HQ  -  mD  x  vQ 


(B-2) 


(B-3) 


Equation  (B-3)  sorves  to  transfer  the  angular  momentum  from  one  reference 
point  0  to  another  reference  point  P. 

If  the  system  is  free  from  external  force,  H'0  and  G  are  constant 
vectors.  Then,  by  Eq.  (B-2),  Hf)  is  a  constant  vector.  In  particular,  if 
ll0  ~  B  and  G  -  0,  it  follows  that  II.,  =  0,  where  P  is  any  point. 


APPENDIX  C 

REMARKS  ON  VECTOR  ANALYSIS 


In  this  Appendix,  n  brief  treatment  of  vector  analysis  is 
presented,  for  more  details,  see  Reference  2,  Chapters  1,  6,  8  and  15. 

In  this  report,  a  letter  with  a  bar  over  it  denotes  a  vector.  Por  example, 

F  stands  for  the  direction  and  magnitude  of  a  force,  although  it  does  not 
signify  the  point  at  which  the  foi'ce  acts.  Ordinarily,  the  point  of  action 
of  a  force  is  designated  by  a  statement.  By  definition,  the  letter  F 
represents  only  the  magnitude  of  vector  F.  Consequently,  F  is  a  non¬ 
negative  number.  A  letter  with  a  caret  over  it  denotes  a  vector  of  unit 

A 

magnitude.  For  example,  F  designates  the  direction  of  force  F . 

The  vector  equation  A  -  B  signifies  that  vectors  A  and  F  have 
the  same  magnitude  and  the  same  direction,  but.  not  necessarily  the  same 
point  of  action.  The  vector  -F  is  defined  to  have  the  same  magnitude  as 
vector  F,  but  the  opposite  direction.  The  vectors  F  and  ~F  need  not  have 
the  same  point  of  action.  For  example,  if  F  denotes  a  force  that  acts  on  a 
body,  the  reaction  of  the  force  is  -F. 

If  k  is  a  positive  number,  the  product  kF  is  defined  to  be  a 
vector  with  the  direction  of  F,  and  with  magnitude  kF.  If  k  is  a  negative 

number,  the  product  kF  is  defined  to  be  a  vector  with  the  direction  of  -F, 
and  with  magnitude  |k|F,  where  |k|  is  the  absolute  value  of  k.  In  view  of 

_  A 

these  definitions,  F  =  FF. 

The  resultant  F  of  two  vectors  F^  and  F^  is  called  the  sum  of  the 

vectors.  The  process  of  obtaining  the  sum  or  resultant  of  two  vectors  by 

the  well-known  parallelogram  construction  is  called  vector  addition. 

Symbolically,  F  =  Fj  +  F2  =*  F2  +  Fj .  It  is  to  he  observed  that  the  relation 

F  =  Fj  +  F2  does  not  imply  that  F  =  +  F2»  In  general,  Fj  +  F2  is  greater 

than  F,  since  the  three  vectors  F,  F^,  F2  form  the  sides  of  a  triangle. 

By  repeated  applications  of  the  parallelogram  construction,  one 

obtains  the  polygon  construction  for  the  sum  of  a  number  of  vectors .  Por 

example,  if  forces  P^ ,  P2,  F^,  act  at  a  point  P  of  a  body,  their 

resultant  F  is  obtained  by  arranging  the  vectors  F^ ,  \:7,  Fj,  F^  in  a  chain, 

7 -  _ 

This  convention  is  not  always  used  in  this  report.  For  example,  M  denotes 

the  moment  of  a  force,  but  M  stands  for  the  mass  of  the  pun.  No  confusion 
should  occur,  since  notations  are  explained. 
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maintaining  the  proper  directions  of  the  vectors.  Then  the  line  segment 
from  the  initial  point  P  to  the  terminal  point  of  the  chain  represents  the 
resul"  .nt  force.  Figure  C-l.  Symbolically,  F  -  F^  +  +  F^.  Vector 

addition  is  commutative  and  associative;  i.e.,  the  order  of  vectors  F^ , 

Fj,  F^  in  the  polygon  is  irrelevant,  and  any  subset  of  the  vectors  in  the 
polygon  may  be  replaced  by  their  resultant.  Symbolically,  this  means  that 
the  vectors  in  the  sum  +  ^  +  ^3  +  ^4  ma^  Permutei^  any  way>  and 
parentheses  may  be  introduced  arbitrarily  in  the  sum. 

Subtraction  of  a  vector  is  defined  to  be  addition  of  the  negative 
vector,  i.e.,  A  -  B  =  A  +  (-B) .  Accordingly,  to  subtract  a  vector  B  from  a 
vector  A,  we  reverse  the  direction  of  B,  and  add  the  reversed  vector  to  A. 

A  direction  can  be  defined  only  with  respect  to  some  reference 
frame.  Consequently,  vector  analysis  cannot  be  entirely  divorced  from 
coordinate  systems.  For  definiteness,  only  right-handed  systems  of  coordi¬ 
nates  are  considered.  ~f  (x,y,z)  is  a  right-handed  system  of  rectangu1 
coordinates,  the  thumb,  the  forefinger,  and  the  middle  finger  of  the  right 
hand  can  be  directed  in  the  positive  senses  along  the  x,  y,  and  z  axes, 
respectively.  Frequently  it  is  convenient  to  designate  the  directions  of 

AAA 

rectangular  coordinate  axes  (x,y,z)  by  three  unit  vectors  (i,j,k)  that 
coincide  in  direction  with  these  axes,  Figure  C-2.  Then  the  position  vector 
R  from  the  origin  to  the  point  P:(x,y,z)  satisfies  the  vector  equation, 

R  =  ix  +  jy  +  kz  (C-l) 

The  orthogonal  projections  (F^jF^jF^)  of  any  vector  F  on  the  axes  (x,y,z) 
satisfy  the  equation 

F  =  iFt  +  jF2  +  kFj  (C-2) 

Addition  and  subtraction  of  vectors,  and  multiplication  of  vectors 
by  scalars  conform  to  the  axioms  of  elementary  algebra.  Consequently, 
cert  operations  with  vectors  can  be  performed  as  in  scalar  algebra.  For 
example,  algebraic  reduction  of  the  vector  equation 

F  -  3 [A  -  Z(B  -  A)  +  5A  -  7B] 
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The  angle  between  a  vector  and  a  cartesian  coordinate  axis  is 
called  the  direction  angle  of  the  vector  with  respect  to  the  coordinate 
axis.  The  three  direction  angles  (a,3,y)  of  a  vector  with  respect  to  three 
rectangular  cartesian  axes  (x,y,z)  determine  the  direction  of  the  vector. 
Figure  C-3.  The  direction  angles  of  a  vector  are  specified  to  lie  in  the 
range  0°  to  180°,  inclusive.  Consequently,  a  direction  angle  is  determined 
uniquely  by  its  cosine.  If  the  cosine  is  negative,  the  angle  is  greater 
than  90°.  The  cosines  of  the  direction  angles  of  a  vector  are  called  the 
direction  cosines  of  the  vector. 


If  (a,8,Y)  are  the  direction  angles  of  a  vector  F,  and  if  F 
denotes  the  magnitude  of  F, 


Fj  =  Fcosa  ,  F2  =  Fcos  B  ,  F^  =  Fcos  y 


(C-3) 


A  direction  in  space  may  be  designated  by  a  unit  vector.  If  F  is  a  unit 
vector,  F  =  1.  Accordingly,  Eq.  (C-3)  shows  that  the  orthogonal  projections 
of  a  unit  vector  on  the  x,  y,  and  z  axes  are  identical  to  the  direction 
cosines  of  the  vector. 

By  trigonometry. 


P2  .  *  F22  4  F32 


(C-4) 


Equations  (C-3)  and  (C-4)  yield 


2  2  2 
cos  a  +  cos  g  +  cos  y  =  1 


(C-5) 


Let  A  and  B  be  two  vectors  whose  projections  on  rectangular 
cartesian  axes  (x,y,z)  are  (Aj^.Aj)  and  (Bj.B^B^).  The  expression 
AjBj  +  ^2^2  +  *s  Ca^ec^  the  scalar  product  (or  dot  product)  of  the 

two  vectors.  This  expression  is  conventionally  denoted  by  A  *  F.  It  is 
seen  by  this  definition  that  A  •  F  =  F  •  A.  Accordingly,  scalar  multiplica¬ 
tion  of  vectors  is  said  to  be  commutative.  The  importance  of  the  scalar 
product  arises  fro r,  a  geometric  identity  that  is  expressed  by  the  equation. 


(C-6) 


A  •  B  =  Aj^Bj  +  A2B2  +  A^Bj  =  ABcosS 

where  A  and  B  are  the  magnitudes  of  vectors  A  and  B,  and  9  is  the  angle 
between  these  vectors. 

Several  special  cases  are  to  be  noted.  If  A  =  F,  Eq.  (C-6)  yields 

2  2  2  2  _ 

A  =  Aj^  +  a2  +  A3  >  which  is  equivalent  to  Eq.  (C-4)  .  If  B  is  a  unit  vec¬ 
tor  (B  =  1) ,  it  is  apparent  from  Eq.  (C-6)  that  A  •  F  is  the  orthogonal 
projection  of  vector  A  on  a  line  with  the  direction  and  sense  of  vector  F. 

If  A  and  B  are  both  unit  vectors,  their  projections  on  the  (x,y,z)  axes  are 
identical  to  their  direction  cosines.  Hence, 

cos  cos  a2  +  cos  cos  B2  +  cos  Yj  cos  y2  =  cos  6  (C-7) 

where  and  (a2>32  »  Y2)  are  the  direction  angles  of  vectors  A  and 

B,  and  9  is  the  angle  between  these  vectors.  If  A  /  0  and  F  /  0,  but 
A  •  B  =0,  vectors  A  and  B  are  perpendicular  to  each  other  ( 9  =  90°). 

Occasionally  an  expression  of  type  A(F  •  C)  arises.  The 
parentheses  may  be  removed;  i.e., 

A(B  ♦  C)  =  A  F  •  C 

There  is  no  ambiguity  in  the  expression  A  B  •  C,  since  no  meaning  is  here 
assigned  to  the  expression  A  F  standing  alone.  Hence,  A  F  •  C  is  a  vector 
with  the  direction  of  vector  A  and  with  magnitude  ABCcos  9  ,  where  0  is  the 
angle  between  vectors  B  and  C. 

Another  expression  that  sometimes  arises  is 

-  AA  -  /S  A  _ 

ii  •  A  +  jj  •  A  +  kk  •  A 

Since  the  (x,y,z)  components  of  A  are  A.  =  i  •  A,  A  =  j  •  A,  and  A_  =  k  •  A, 

1  Z  o 

this  reduces  to 


iAj  +  jA2  +  kAj  =  A 
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Consequently, 


ii  •  A  +  jj  •  A  +  kk  •  A  =  A  (C-8) 

This  relation  is  an  identity. 

By  cartesian  expansion,  it  is  easily  seen  that 

A  •  (F  +  C)  =  A  •  F  +  A  •  C 

Since  also  A  •  B  =  B  •  A,  the  scalar  product  conforms  to  the  rules  of 
elementary  algebra. 

From  two  given  vectors,  A  and  B,  a  third  vector  C  may  be  derived 
by  the  definition, 

C1  =  A2B3  ~  A3B2  »  C2  =  A3B1  “  ^l^  »  C3  =  AiB2  "  A2B1  (c“9) 
This  may  be  expressed  concisely  in  determinant  notation: 


A 

i 

C  =  Ax 

A 

j 

A2 

k 

A3  =  A  x  B 

(C-10 

B1 

B2 

B3 

The  vector  C  is  called  the  vector  product  or  cross  product  of  A  and  B.  It 
may  be  shown  by  geometry  that  the  magnitude  of  C  is 

C  =  ABsinQ  (C-ll) 

where  0  is  the  angle  between  vectors  A  and  B.  It  follows  from  Eq.  (C-9)  that 
C,A=C,B=0.  Consequently,  vector  C  is  perpendicular  to  both  of  the 
vectors  A  and  B.  It  can  be  shown  that,  if  the  coordinates  (x,y,z)  are 
right-handed,  the  sense  of  vector  C  is  that  in  which  a  right-hand  screw 
advances  when  turned  from  A  to  B. 

The  vector  product  is  not  commutative.  Since  a  permutation  of 
two  rows  in  a  determinant  changes  the  sign  of  the  determinant,  Eq.  (C- 10) 


shows  that  B  x  A  =  -A  x  B.  In  spite  of  this  anomalous  behavior,  the  vector 
product  has  other  properties  in  common  with  ordinary  multiplication.  In 
particular, 

R  x  (A  +  B)  sRjA  +  RxB,  (A+B)  xR-AxR  +  BxR 

(C- 12) 

Hence , 

(A+B)  x(C  +  D)  =  AxC  +  AxD  +  BxC  +  BxD 

It  is  seen  from  Eq.  (C-11)  that  the  vector  product  of  two  parallel  vectors 
is  zero,  since  0=0.  Hence, 

A  A  A  A  A  A 

ixi  =  jxj=kxk  =  0 
Also,  by  Eq.  (C-10) 

AAA 

i  j  k 

A  A  A 

l  x  j  =  1  0  0  =  k 

0  10 

Similarly,  k  x  i  =  j  and  j  x  k  =  i .  Evidently, 

A  x  B  =  (iAj  +  jA2  +  £a3)  x  (lBj  +  jB2  +  kB^) 

Algebraic  expansion  of  the  right  side  of  this  equation  leads  back  to 
Eq.  (C-9) . 

The  expression  A  •  (B  x  C)  is  called  the  scalar  triple  product . 

It  may  be  written  without  parentheses  as  A  •  B  x  C,  since  (A  •  B)  x  C  has 
no  meaning.  The  expression  A  •  B  x  C  is  a  scalar.  Cartesian  expansion 
yields  the  determinant  form, 


(C-13) 


A1  A2  A3 

A  •  B  x  C  =  Bj  B2  Bj 

C1  C2  S 


Since  a.  transposition  of  two  rows  in  a  determinant  merely  changes  the  sign  of 
the  determinant, 


A*  B  xC  =  B  •  C  x  A  5  C  •  AxB 


CC-14) 


The  absolute  value  of  A  •  B  x  C  represents  the  volume  of  the  parallelepiped 
having  concurrent  edges  represented  by  A,  B,  and  C. 

The  vector  triple  product  is  A  x  (B  x  C) .  The  parentheses  are 
essential  in  this  expression.  By  cartesian  expansion,  the  following  identity 


can  be  verified: 


A  x  (B  x  C)  =  B  A  •  C  •  C  A  •  B 


(C-15) 


This  may  be  memorized  as  the  "Back-Cab"  formula. 

The  vector  product  is  useful  for  representing  moments  of  forces . 
Let  a  be  the  perpendicular  distance  from  a  given  point  0  to  the  line  of 
action  of  a  given  force  F,  Figure  C-4 .  The  moment  of  force  F  about  point  0 
is  defined  to  be  a  vector  M  with  magnitude  Fa.  The  vector  M  is  defined  to  be 
perpendicular  to  the  plane  determined  by  the  force  F  and  point  0.  The  sense 
of  vector  M  is  defined  by  the  right-hand-screw  rule;  i,e.,  vector  M  points  in 
the  direction  that  a  right-hand  screw  would  advance  if  force  F  should  cause 
it  to  turn  about  an  axis  through  point  0.  For  the  case  illustrated  by 
Figure  C-4,  vector  M  is  directed  toward  the  reader,  perpendicular  to  the 
plane  of  the  paper.  Let  the  vector  OP  (Figure  C-4)  be  denoted  by  r.  Then 
the  conditions  of  the  preceding  definition  are  fulfilled  by  the  vector 


equation. 


M  =  r  x  F 


(C-16) 
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The  moment  of  a  force  about  an  axis  is  a  scalar.  If  n  is  a  unit 
vector  in  the  direction  of  axis  L,  the  moment  of  force  F  about  axis  L  is 
•  M,  whei*e  M  is  the  moment  of  force  F  about  any  point  0  on  axis  L. 
Consequently,  by  Eq.  (C-16)  , 

ML=F*Axr=r-Fxn=fi*rxF  (C-17) 

An  infinitesimal  increment  dR  of  a  vector  R  need  not  be  collinear 
with  the  vector  R,  Figure  C-5.  Consequently,  in  general,  the  vector  R  +  dR 
differs  from  the  vector  R,  not  only  in  magnitude,  but  also  in  direction. 

It  would  be  misleading  to  denote  the  magnitude  of  vector  dR  by  dR,  since  dR 
denotes  the  increment  of  the  scalar  R.  Accordingly,  the  magnitude  of  dR  is 
denoted  by  | dR |  ,  or  by  another  symbol,  such  as  ds .  The  magnitude  of  the 
vector  R  +  dR  is  R  +  dR.  Figure  C-S  shows  that  dR  £  |dR| .  If  the  vector  R 
is  a  function  of  a  scalar  t  (where  t  may  or  may  not  denote  time) ,  dR/dt  is 
defined  to  be  a  vector  in  the  direction  of  dR  with  magnitude  ds/dt,  where 

ds  =  | dR | .  If  R  is  the  position  vector  of  a  particle,  and  if  t  denotes 

—  —  —  2—  2 
time,  v  =  dR/dt  is  the  velocity  of  the  particle  and  dv/dt  =  d  R/dt  is  the 

acceleration  of  the  particle.  Vectors  obey  the  same  rules  of  differentiation 

as  scalars.  This  may  be  shown  by  the  delta  method  that  is  used  for  deriving 

differentiation  formulas  in  scalar  calculus.  For  example,  if  Q  =  uR,  where 

u  is  a  scalar  function  of  t  and  R  is  a  vector  function  of  t, 

dO  —  •  —  -- 

=  Q  =  uR  +  uR 

in  which  the  dot  denotes  the  derivative  with  respect  to  t.  Likewise, 

^(A  •  B)  =  A  •  B  +  A  •  B 

and  ~(K  xB)  =  AxB+AxB 

The  angular  velocity  w  of  a  rigid  body  is  a  vector  quantity.  For, 
let  u)  represent  the  angular  velocity  of  a  rigid  body  whose  motion  at  the 
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instant  under  consideration  is  a  rotation  about  an  axis  L.  Let  r  be  the 
vector  from  any  point  0  on  axis  L  to  a  particle  P  of  the  body  (Figure  C-6) . 
The  distance  of  particle  P  from  the  axis  of  rotation  is  rsinG  ,  where  9  is 
the  angle  between  the  vectors  r  and  a>.  Hence,  the  speed  of  particle  P  is 
rt»)sin  9  .  The  velocity  v  of  particle  P  is  perpendicular  to  the  plane  of 
vectors  r  and  w,  and  its  sense  is  determined  by  the  right-hand-screw  rule 
(if  right-handed  coordinates  are  used) .  Therefore, 

v  =  oJ  x  r  (C-18) 
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ATTN:  Dr.  M.  Beckstead 
Provo,  UT  84601 

1  University  of  California 
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P.O.  Box  808 
Livermore,  CA  94550 
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Dept,  of  Mechanical  Engineering 
ATTN :  Dr.  P.  C.  Chou 
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DRXSY-LM,  J.C.C.  Fine 
Dir,  USAHEL 

ATTN:  DRX1IE ,  Dr.  J.D.  Weisz 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1.  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  pi*epared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


FOLD  HERE 


Director 

US  Army  Ballistic  Research  Laboratory 

ATTN:  DRDAR-BIA-S 

Aberdeen  Proving  Ground,  MD  21005 
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Director 

US  Army  Ballistic  Research  Laboratory 

ATTN:  DRDAR-  BLA-S 

Aberdeen  Proving  Ground,  Ml)  2100$ 


FOLD  HERB 


